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B.  TECHNICAL  SECTION 

I .  Abstract 

x  The  overall  objective  of  this  proposal  is  to  investi¬ 
gate  the  robustness  to  departures  from  independence  of  methods 
currently  in  use  in  reliability  studies  when  competing  failure 
modes  or  competing  causes  of  failure  associated  with  a  single 
mode  are  present  in  a  series  system.  The  first  specific  aim 
is  to  examine  the  error  one  makes  in  modeling  a  series  system 
by  a  model  which  assumes  statistically  independent  component 
lifetimes  when  in  fact  the  component  lifetimes  follow  some 
multivariate  distribution.  The  second  specific  aim  is  to  assess 
the  effects  of  the  independence  assumption  on  the  error  in 
estimating  component  parameters  from  life  tests  on  series 
systems.  In  both  cases,  estimates  of  such  errors  will  be  deter¬ 
mined  via  mathematical  analysis  and  computer  simulations  for 
several  prominent  multivariate  distributions.  A  graphical  dis¬ 
play  of  the  errors  for  representative  distributions  will  be 
made  available  to  researchers  who  wish  to  assess  the  possible 
erroneous  assumption  of  independent  competing  risks.  A  third 
aim  is  to  tighten  the  bounds  on  estimates  of  component  relia¬ 
bility  when  the  risks  belong  to  a  general  dependence  class  of 
distributions  (for  example ,  positive  quadrant  dependence,  posi¬ 
tive  regression  dependence,  etc.)v  Major  decisions  involving 
reliability  studies,  based  on  competing  risk  methodology,  have 
been  made  in  the  past  and  will  continue  to  be  made  in  the  future 
This  study  will  provide  the  user  of  such  techniques  with  a 
clearer  understanding  of  the  robustness  of  the  analyses  to  de¬ 
partures  from  independent  risks,  an  assumption  commonly  made 
by  the  methods  currently  in  use.  \ 
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II.  Specific  Objectives 

The  overall  objective  is  to  investigate  the  robust¬ 
ness  to  departures  from  independence  to  methods  currently  in 
use  in  reliability  studies  when  competing  failure  modes  or  com¬ 
peting  causes  of  failure  associated  with  a  single  mode  are 
present  in  a  series  system.  We  shall  also  refer  to  such  com¬ 
petitive  events  as  competing  risks.  The  approach  will  be  through 
the  investigation  of  certain  aspects  of  specific  parametric  multi 
variate  distributions  or  by  classes  of  distributions  which  are 
appropriate  in  reliability  analyses  when  there  are  competing 
risks  present. 

The  specific  objectives  are: 

» 

1)  to  assess  the  error  incurred  in  modeling  system 
life  in  a  series  system  assumed  to  have  indepen¬ 
dent  component  lifetimes  when  in  fact  the  com¬ 
ponent  lifetimes  are  dependent. 

2)  to  assess  the  error  in  estimating  component  param¬ 
eters  (i.e. ,  component  reliability,  mean  com¬ 
ponent  life,  etc.)  in  a  series  system  employing 
either  parametric  or  nonparametric  models  which 
assume  independent  component  failure  times  when 

in  fact  the  lifetimes  are  dependent  and  follow 
some  plausible  multivariate  distribution.*  . 

3)  to  derive  bounds  on  component  reliability  when 
the  failure  modes  are  dependent  and  fall  in  a 
particular  dependence  class  (e.g.,  positive  quad¬ 
rant  dependence,  positive  regression  dependence, 
etc. ) . 

4)  to  develop  tests  of  independence,  based  on  data 
collected  from  series  systems,  by  making  some 
restrictive  assumption  about  the  structure  of  the 
systems. ** 

*  A  plausible  parametric  multivariate  distribution  will  be 
one  that  satisfies  one  of  the  following  conditions: 

i)  the  distribution  of  the  minimum  of  the  component 
failure  times  closely  approximates  widely  accept¬ 
ed  families  of  system  life  distributions. 

or  ii)  the  marginal  distributions  closely  approximate 
the  distributions  of  component  failure  times  in 
the  absence  of  other  failure  modes. 


**This  objective  has  been  added  to  the  original  objectives  be¬ 
cause  it  answers  a  natural  question  raised  by  our  preliminary 
investigation. 
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XXX.  *  Introduction  to  Problem  and  Significance  of  Study 

•  .  • 

•  •  •  - 

Alvin  Weinberg  (1978)  in  an  editorial  comment  in  the ' 
published  proceedings’  of  a  workshop  on  Environmental  Biologi¬ 
cal  Hazards  and  Competing  Risks  noted  that  "the  question  of 
competing  risks  will  not  quiet lyt  go  away:  corrections  for  com¬ 
peting  risks  should  be  applied  routinely  to  data.*  The  problem 
of  competing .  risks  commonly  arises  in  .a  wide  range  of  experi¬ 
mental  situations.  Although  we  shall  confine  our  attention 
in  the  following  discussion  to  those  situations  involving 
series  systems  in  which  competing  failure  modes  or  competing 
causes  of  failure' associated  with  a  single  mode  are  present, 
it  is  certainly  true  that  we  might  just-  as  easily  speak  of 
clinical  trials,  animal  experiments,-  or  other  medical  and  bio¬ 
logical  a  todies  where  competing  events  interrupt  our  study  of 
the  main  event  of  interest  (cf  •  Lagokos  (1979))  • 

Consider  electronic  or  mechanical  systems,  such  as 

i *~t  transmission  equipment,-  computers,  aircraft,  missiles 
e«d  other  weaponry  consisting  of . several  components  in  series. 
Usually  each  component  will'  have  a  random  life  length  and  the 
life  of  the  entire  system  will  end  with  the  failure  of  the 
shortest  lived  component.  We  will  examine  two  situations  more 
closely  in  which  competing .  risks  play  a  vital  role. 

First,  suppose  we  are  ‘  attempting  to  evaluate  system  life 
from  knowledge  of  the  individual  component  lifetimes.  Such 
an  evaluation  will  utilize  either  an  analysis  involving  math¬ 
ematical  statistics  or  a  computer  simulation. ..  At  a  recent 
conference- 'an  Modeling  and  Simulation,  McLean  (1981)  presented 
a  scheme  to  simulate  the  life  of  a 'missile  Which  consisted  of 
many  major*  components  in  series.  -  The  failure  distribution  asso¬ 
ciated  with  each- component  was  assumed  to  be  known  (usually 

*0  or  WeibulT. )  To  arrive  at  the  system  failure  dis¬ 
tribution,  the  components  were-  assumed  to  act  independently  of 
each  other.  Realistically,  this  may  or  may  not  be  the  case. 

Xf  the  component  lifetimes  were  dependent  for  any  reason,  the 
computed  system  failure  distribution  (as  well  as  its  subsequent 
parameters  such  'as  system  mean  life  and  system  reliability  for 
a  specified  time)  would  only  crudely  approximate  the  true 
distribution.  The  first  specific  aim  of  this  proposal  is'  to 
ascertain  the  ^rrc***  incurred  in  modeling  system  life  in  a 
series  system  assumed  to  have  independent  component  lifetimes 
(i.e.,  risks)  when,  in  fact,  the  risks  are  dependent. 

Second,  suppose  we  wish  to  evaluate  same  aspect  of  the 
distribution  of  a  particular  failure  mode  based  on  a  typical 
life  test  of  a  series  system.  The  response  of  interest  is  the 
time  until  failure  of  a  particular,  mode  of  interest.  Frequently 
this  response  will  hot  be  observable  due  to  the  occurrence  of 
some  other  event  which  precludes  failure  associated  with  the 
mode  of  interest.  We  shall  term  such  competing  events  which 
interrupt  our  study  of  the  failure  modes  of  interest  as 

competing  risks. 
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Competing  risks  arise  in  such  reliability  studies  when 

1)  the  study  is  terminated  due  to  a  lack  o£  funds  or  the 
pr e-de termined  period  of  observation  has  expired 
(Type  Z  censoring) . 

2)  the  study  is  terminated  due  to  a  pre-de termined  number 
of  failures  of  the  particular  failure  mode  of  interest 
being  observed  (Type  XZ  censoring) . 

3)  some  systems  fail  because  components  other  than  the 
one  of  interest  malfunction. 

4)  the  component  of'  interest  fails  from  some  cause  other 
than  the  one  of  interest. 

*  • 

Xn  all  four  situations,  one  may  think  of  the  main  event  of 
interest  as  being  censored,  i.e.,  not  fully  observable.  In  the 
first  two  situations,  the  time  to  occurrence  of  the  event  of 
interest  should  be  independent  of  the  censoring  mechanism,  in 
such  instances,  the  methodology  for  estimating  relevant  reliabili¬ 
ty  probabilities  has  received  considerable  attention  (cf.  David 
and  Moeschberger  (1978)',  Kalbfeish  and  Prentice  (1980) ,  Eland t- 
johnson  and  Johnson  (1980) ,  Mann,  Schafer,  Singpurwalla  (1974) 
and  Barlow  and  Proschan  (1975)  for  references  and  discussion)' - 
In  the  third  situation,  the  time  to  failure  of  the  component  of 
interest  may  or  may  not  be  independent  of  the  failure  times  of 
-other  components  in  the  system.  For  example,  there  may  be 
common  environmental  factors  such  .as  extreme  temperature  which 
may  affect  the  lifetime  of  several  components.  Thus  the  question 
of  dependent'  coupe  ting  risks  is  raised.  A  -similar  observation 
may  be  made  with  respect  to  the  fourth  situation,  viz.,  failure 
times  associated  with  different  failure  modes  of  a  single  com¬ 
ponent  may  be  dependent.  For  a  very  special  type  of  dependence, 
the  models  discussed  by  Marshall-Olkin  (1967),  Langberg,  Proschan 
and  Quinsy  (1978)  ,  and  Langberg,  Proschan,  and  Quinsy  (1981) 
allow  one  to  convert  dependent  models  into  independent  ones. 

Zf  no  assumptions  whatever  are  made  about  the  type  of  * 
dependence  between  the  distribution  of  potential  failure  times, 
there  appears  to  be  little  hope  of  estimating  relevant  component 
parameters.  In  some  situations,  one  may  be  appreciably  misled 
(cf.  Tsiatis  (1975),  Peterson  (1976)).  However,  as  Easterling 
(1980)  so  clearly  points  out  in  his  review  of  Bimbaum's  (1979) 
monograph 

"there  seems  to  be  a  need  for  some  robustness 
studies.  How  far  might  one  be  off,  quantita¬ 
tively,  if  his  analysis  is  based  on  incorrect 
assumptions? " 

The  second  specific  aim  will  address  this  important 
issue.  First  if  a  specific  parametric  model  which  assumes 
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independent  risks  has  been  used  in  the  analysis,  it  would 
be  of  interest  to  know  how  the  error  in  estimation:  is 
affected  by  this  assumption  of  independence.  That  is,  if 
independent  specific  parametric  distributions  are  assumed 
for  the  failure  times  associated  with  different  failure 
modes  when  we  really  should  use  a  bivariate  (or  multivariate) 
distribution,  then  what  is  the  .magnitude  of  the.  error  in 
estimating -component  parameters?  Secondly,  one  may  wish  to 
allow  for  a  less  stringent  type  of  model  assumption,  and  ask 
the  same  question  with  regard  to  the  estimation  error.  That  ' 
is,  if  a  nonparametric  analysis  is  performed,  assuming  in¬ 
dependent  risks,  when  some  types  of  dependencies  may  be 
present,  then  what  is  the  magnitude  of  the  estimation  error? 

The  third  specific  aim  will  attempt  to  obtain  bounds  on 
the  component  reliability  when  the  failure  times  belong  to 
a  broad  dependence  class  (e.g. ,  association,  positive  quadrant . 
dependence,  positive  regression  dependence,,  etc.}*  More 
details  will  be  presented  in  the  methods  section.' 

In  summary,  competing  risk  analyses  have  been  performed 
in  the  past-  will  continue  to  be  performed  in  the  future. 
This  study  will  provide  the  user  of  such  techniques  with  a 
clearer  understanding  of  the  robustness  to  departures  from 
independent  risks,  an  assumption  which  most  of  the  methods 
currently  in  use  assume. 
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IV.  Progress  Report  on  Third  Year's  Work 

A  summary  of  the  first  and  second  year's  work  is  reported  in  the  annual  reports 
dated  October  26th  1983  and  October  26th  1984  respectively.  We  believe  that  during  the 
past  three  years  we  have  made  substantial  progress  in  dealing  with  the  objectives  as 
outlined  on  page  4.  In  addition  to  the  papers  and  articles  referred  to  in  the  first  two 
annual  reports,  we  would  like  to  mention  some  of  the  more  recent  work. 

First,  the  recently  published  paper  which  investigates  the  problem  of  improving  the 
product-limit  estimator  of  Kaplan  and  Meier  (1958)  when  there  is  extreme  independent 
right  censoring  is  presented  in  Appendix  A.  This  paper  looks  at  several  techniques  for 
completing  the  product  limit  estimator  by  estimating  the  tail  probability  of  the  survival 
curve  beyond  the  largest  observed  death  time.  Two  methods  are  found  to  work  well  for  a 
variety  of  underlying  distributions.  The  first  method  replaces  those  censored 
observations  larger  than  the  biggest  death  time  by  the  expected  order  statistics, 
conditional  on  the  largest  death,  computed  from  a  Weibull  distribution.  The  WeibuU  is. 
chosen  since  it  is  known  to  be  a  reasonable  model  for  survival  in  many  situations. 
Parameters  of  the  model  are  estimated  in  several  ways,  but  the  method  of  maximum 
likelihood  seems  to  provide  the  best  results.  The  second  method  replaces  the  constant 
value  of  the  product  limit  estimator  beyond  the  last  death  time  by  the  tail  of  a  Weibull 
survival  function.  Again  parameters  are  estimated  by  a  variety  of  methods  with  the 
maximum  likelihood  estimators  performing  the  best. 

Second,  a  paper  which  obtains  bounds  on  the  component  reliability,  based  on  data 
from  a  series  system,  for  the  Oakes  (I  982)  model  has  been  revised.  Since  this  model  has 
the  same  dependence  structure  as  the  random  effects  model  with  w  having  a  gamma 
distribution,  these  bounds  are  good  for  a  general  class  of  distributions.  The  bounds, 
which  are  determined  by  specification  of  a  range  of  coefficients  of  concordance,  are 
found  by  solving  a  differential  equation  in  the  observable  system  reliability  and  crude  life 
on  one  hand  and  the  unobservable  component  survival  function  on  the  other  hand.  This 
revision  is  reproduced  in  Appendix  B. 
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Third,  we  have  submitted  an  overview  paper  for  publication  which  summarizes 
some  of  the  work  performed  during  the  past  three  years.  The  results  of  this  paper  were 
presented  in  an  invited  talk  to  the  Eastern  North  American  Region  of  the  Biometrics 
Society  at  Raleigh,  North  Carolina  in  the  Spring  of  1985.  (See  Appendix  C  for  a  copy  of 
this  paper.) 

Finally,  a  paper  has  been  developed  which  discusses  some  general  properties  of  a 
random  environmental  stress  model.  Estimation  of  parameters  under  the  Gamma  stress 
model  is  considered,  and  a  new  estimator  based  on  the  scaled  total  time  on  test 
transform  is  presented  (See  Appendix  D).  Part  of  the  results  in  this  paper  were  presented 
at  the  International  Statistical  Institute  meeting  in  Amsterdam  in  August,  1985.  A  copy 
of  that  contributed  paper  is  found  in  Appendix  E. 

V.  Methods 

We  refer  to  pages  8-52  of  the  original  proposal  for  a  discussion  of  the  general 
methodology. 
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A  Comparison  of  Several  Methods  of  Estimating  the  Survival 
Function  when  There  Is  Extreme  Right  Censoring 

M.  L.  Moeschberger1  and  John  P.  Klein2 

Departments  of  Preventive  Medicine1  and  Statistics2,  The  Ohio  State  University, 
320  West  10th  Avenue,  Columbus,  Ohio  43210,  U.S.A. 


Summary 

When  there  is  extreme  censoring  on  the  right  the  Kaplan-Meier  product-limit  estimator  is  known  to 
be  a  biased  estimator  of  the  survival  function.  Several  modifications  of  the  Kaplan-Meier  estimator 
are  examined  and  compared  with  respect  to  bias  and  mean  squared  error. 


1.  Introduction 

In  human  and  animal  survival  studies,  as  well  as  in  life-testing  experiments  in  the  physical 
sciences,  one  method  of  estimating  the  underlying  survival  distribution  (or  the  reliability 
of  a  piece  of  equipment)  which  has  received  widespread  attention  is  the  Kaplan-Meier 
product-limit  estimator  (Kaplan  and  Meier,  19S8). 

For  the  situation  in  which  the  longest  time  an  individual  is  in  a  study  (or  on  test)  is  not 
a  failure  time,  but  rather  a  censored  observation,  it  is  well  known  that  there  are  many 
complex  problems  associated  with  any  statistical  analysis  (Lagakos,  1979).  In  particular, 
the  Kaplan-Meier  product-limit  estimator  is  biased  on  the  low  side  (Gross  and  Clark, 
1975).  In  the  case  of  many  censored  observations  larger  than  the  largest  observed  failure 
time,  this  bias  tends  to  be  worse.  Estimated  mean  survival  time  and  selected  percentiles,  as 
well  as  other  quantities  dependent  on  knowledge  of  the  tail  of  the  survival  function,  will 
also  exhibit  such  biases. 

A  practical  situation  which  motivates  this  study  is  a  large-scale  animal  experiment 
conducted  at  the  National  Center  for  Toxicological  Research  (NCTR),  in  which  mice  were 
fed  a  particular  dose  of  a  carcinogen.  The  goal  of  the  experiment  was  to  assess  the  effects 
of  the  carcinogen  on  survival  and  on  age-specific  tumor  incidence.  Toward  this  end,  mice 
were  randomly  divided  into  three  groups  and  followed  until  death  or  until  a  prespecified 
group  censoring  time  (280,  420,  or  560  days)  was  reached,  at  which  time  all  those  still  alive 
in  a  given  group  were  sacrificed.  Often  there  were  many  surviving  mice  in  all  three  groups 
at  the  sacrifice  times. 

In  general,  we  consider  an  experiment  in  which  n  individuals  are  under  study  and 
censoring  is  permitted.  Let  ;<n,  ■  •  • ,  t{m)  denote  the  m  ordered  failure  times  of  those  m 
individuals  whose  failure  times  are  actually  observed  (4i>  <  ■  <  The  remaining 

n  —  m  individuals  have  been  censored  at  various  points  in  time.  It  will  be  useful  to  introduce 
the  notation  5,  to  denote  the  number  of  survivors  just  prior  to  time  UjV,  that  is,  St  is  the 
number  of  individuals  still  under  observation  at  time  r0),  including  the  one  that  died  at 
tij).  Then  the  Kaplan-Meier  product-limit  estimator  (assuming  no  ties  among  the  r(J))  of 


Key  words:  Adjusted  Kaplan-Meier  suro  vaJ  estimation;  Bias  of  survival  function;  Life-testing;  Right 
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the  underiying  survival  function.  P(t)  =■  PriT >  t),  is 

l  for 

P{t)  -  !  n  (■&  "  1  VS,  for 

fm 1 

0  for 


t  <  f<n 

hi >  <  /  <  hj+i) 

t  * 


(1) 


for;  *  l, ....  m,  where  » te  if  the  longest  time  an  individual  is  on  study  is  a  censoring 
time  or  r<m*»  ■  «  if  the  longest  time  an  individual  is  on  study  is  a  death. 

This  paper  first  proposes,  in  §2,  some  methods  of  “completing”  the  Kaplan-Meier 
estimator  of  the  survival  function  by  (i)  replacing  those  censored  observations  that  are 
larger  than  the  last  observed  failure  time  by  their  expected  order  statistics;  (ii)  using  a 
Weibull  distribution  to  estimate  the  tail  probability  P(t),  for  t  >  tc\  and  (iii)  employing  a 
method  suggested  by  Brown,  Hollander,  and  Korwar  (BHK)  (1974).  The  second  purpose 
is  to  demonstrate  the  magnitude  of  the  bias  and  mean  squared  error  (MSE)  of  the  Kaplan- 
Meier  estimator  and  to  compare  all  methods  of  “completing”  P{t)  in  the  context  of  the 
aforementioned  mouse  study,  utilizing  simulated  lifetimes  from  exponential,  Weibull, 
lognormal,  and  bathtub-shaped  hazard  function  distributions.  These  results  are  presented 
in  §3. 


2.  Completion  of  Kaplan-Meier  Product-Limit  Estimator 


2. 1  Expected  Order  Statistics 

One  method  of  attempting  to  “complete”  P{t),  t  >  tc,  would  be  to  “estimate”  the  failure 
times  for  those  censored  observations  that  are  larger  than  the  longest  observed  lifetime.  Let 
nc  be  the  number  of  censored  observations  larger  than  t(m).  A  theorem  regarding  the 
conditional  distributions  of  order  statistics  states  that  for  a  random  sample  of  size  n  from 
a  continuous  parent,  the  conditional  distribution  of  Tiu),  given  T,  «— ne)  m  hi r— rtf)*  u  >  n  nc, 
is  just  the  distribution  of  the  (u  -  n  +  njth  order  statistic  in  a  sample  of  size  drawn  from 
the  parent  distribution  truncated  on  the  left  at  t  -  (see  David,  1981,  p.  20). 

For  computational  purposes,  take  tc  as  an  estimate  of  the  (n  -  «r)th  order  statistic.  Then 
find  the  expected  value  of  the  n(  order  statistics  from  the  parent  distribution  truncated  on 
the  left  at  f<.  Since  the  Weibull  distribution  with  survival  function  P{t)  »  exp (-tk/9)  has 
been  widely  accepted  as  providing  a  satisfactory  fit  for  lifetime  data,  it  seems  reasonable  to 
employ  the  results  of  Weibull  distribution  theory  to  complete  P{t),  t  >  tc.  (It  should  be 
noted  that  any  distribution  which  is  reasonable  for  the  specific  situation  may  be  used.)  The 
expected  values  of  Weibull  order  statistics  up  to  sample  size  40  for  location  parameter 
equal  to  1  and  shape  parameter  equal  to  .5  (0.5)4<  1  )8  may  be  found  in  Harter  ( 1 969).  For 
larger  sample  sizes,  he  states  a  recurrence  relation  which  may  be  used. 

To  compute  expected  values  of  the  n<  order  statistics  in  question,  values  for  ir  and  B  must 
be  chosen.  One  approach  is  to  use  the  maximum  likelihood  estimators,  £  and  8,  computed 
by  using  all  observations  to  estimate  k  and  9.  A  second  approach,  due  to  White  (1969), 
uses  least  squares  estimates  of  k  and  9  obtained  by  fitting  the  model 

ln^^J-d/ZcjlnS-hd/^lnf/fd,;,)]  (2) 

to  the  hj)\  where  H(hj))  is  the  estimated  cumulative  hazard  rate  at  hi)  obtained  from  the 
Kaplan-Meier  estimator.  In  our  Monte  Carlo  study,  we  found  the  maximum  likelihood 
estimators  performed  better  than  the  least  squares  estimators  in  all  cases.  Consequently, 
the  method  of  least  squares  will  be  dropped  from  future  discussion  in  this  paper. 

The  survival  function  for  a  Weibull  random  variable,  truncated  on  the  left  at  tc,  is 

PAD  -  exp(-(f*  -  tkc)/9],  t  >  tc.  (3) 


v-'.vjy 

r ji  r  *  r 
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So,  by  the  theorem  on  order  statistics  stated  at  the  beginning  of  this  section,  the  conditional 
distribution  of  Tfu),  given  (u  ■  n  -  n<  +  1, . . . ,  n)  will  be  approximated 

by  the  (u  -  n  +  n,)th  order  statistic  in  a  sample  of  nc  drawn  from  (3).  For  simplicity,  let 
j  »  u  -  n  +  nt,  so  that  j  -  1, . . . ,  Now  the  expected  value  of  the  yth  order  statistic  from 
(3)  is 

E £  t[Pr<t)rl[PAt)],,'-'+'(ktk-'/8)  dt 

-  /)  X"  ( +  dy  (4) 

where  P{  y)  ■  exp(-y*/tf ),  y  ■  (f*  —  /J)17*  »  0  and  r;:„f  is  the /th  order  statistic  in  a  sample 
of  size  nc.  Equation  (4)  can  also  be  written  as 

E (Tj:„c)  -  ne(*~  ^  j  £  +  tf)'A[F(z)rl[^)]','~;"'^~'  *  (5) 

where  P{z)  «  expC-r*),  z  -  ( y/0)lA  »  0.  Now  E(Tj:„e)  may  be  crudely  estimated  by 

(6) 

where  E(Z/;Br)  is  the  expected  value  of  the  yth  order  statistic  from  a  sample  of  size  nc 
determined  from  Harter’s  (1969)  tables  or  recurrence  relation,  and  8  and  k  are  maximum 
likelihood  estimators  of  8  and  k,  respectively. 

These  nc  estimated  expected  order  statistics  may  then  be  treated  as  “observed”  lifetimes  - 
in  adjusting  (or  “completing”)  the  estimated  survival  function  computed  in  ( 1).  The  area 
under  the  estimated  survival  function  up  to  tc  remains  unchanged.  The  area  under  the 
extended  estimated  survival  function  based  on  the  estimated  expected  order  statistics  is 
then  added  to  the  initial  area  to  obtain  a  more  precise  estimate  of  P(t)  [estimated  order 
statistic  (EOS)  extension]. 

2.2  W exbull  Maximum  Likelihood  Techniques 

A  straightforward  approach  to  completing  P(t)  is  to  set 

P{t)  ■  exp(-f*/9)  for  t  >  tc.  (7) 

Estimates  of  k  and  9  based  on  all  observations  can  be  obtained  by  either  the  maximum 
likelihood  (WTAIL)  or  the  least  squares  method.  However,  our  study  found  the  completion 
using  maximum  likelihood  estimators  was  always  better  in  terms  of  bias  and  mean  squared 
error. 

One  suggestion  for  ostensibly  improving  this  estimator  would  be  to  “tie”  the  estimated 
tail  to  the  product-limit  estimator  at  rc.  Two  methods  were  attempted  to  accomplish  this 
goal  First,  the  likelihood  was  maximized  with  respect  to  k  and  9  subject  to  the  constraint 
that  exp (-r£/0)  -  P(te).  This  method  will  be  referred  to  as  the  restricted  MLE  tail  probability 
estimate  (RWTAIL  extension).  Second,  a  scale-shift  was  performed  on  the  tail  probability 
in  (7)  to  tie  it  to  the  product-limit  estimator.  This  method  led  to  higher  biases  and  mean 
squared  errors  of  the  survival  function  and  will  be  dropped  from  further  discussion  in  this 
paper. 

2.3  BHK-Type  Methods 

The  Brown- Hollander- Korwar  completion  of  the  product-limit  estimator  sets 
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where  9*  satisfies  P(te)  •  exp (-te/Om).  In  the  BHK  spirit  we  tried  to  complete  P(t)  by  a 
Weibuil  function  which  used  estimates  of  k  and  9,  k*  and  9*,  that  satisfied  the  following 
two  relations: 

Pitm)  -  exp(-<&/n 

and 

Pikm- d)  -  exp (~tfc-{)/6m). 

The  latter  method  also  led  to  consistently  poor  performance  and  the  results  will  not  be 
presented. 

Tahiti 


Bias/ 100  {and  MSB/ 100*)  for  estimating  mean  survival  time  for  various  methods  of  completion 
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.208 

299 

(7.886)’ 

(4.498) 

(1.922)* 

(1344) 

(3292) 

600 

25.5 

-3.625- 

-1704 

-.187* 

.344 

.479 

(13.179)- 

(7_522) 

(3.025)* 

(4.275) 

(6.031) 

400 

24.6 

-.991- 

-.047 

-.046 

.016* 

.0379 

(1.011)- 

(.215)* 

(157) 

(175) 

(-343) 

km  1 

500 

316 

-1.632- 

-.049 

-.047* 

.073 

.116 

(1696)* 

(.416)* 

(.535) 

(.508) 

(.705) 

600 

39.3 

-1359- 

.022* 

.034 

.140 

214 

(5.592J* 

(.596)* 

(.987) 

(1.023) 

(1.353) 

400 

7.5 

-.036 

.136* 

-.005 

.003* 

.004 

(.012)* 

(.053)- 

(.013) 

(.014) 

(.014) 

* 

« 

500 

34.6 

-.314 

1.507* 

-.020 

.014* 

.019 

(.109) 

(2830)- 

(.036)* 

(.041) 

(.044) 

600 

59.9 

-.903 

5.982- 

.144 

.028* 

.039 

(.822) 

(41.430)- 

(4.168) 

(.147)* 

(.157) 

Lofnormal 

400 

20.6 

-.868* 

-.178* 

-.544 

-.586 

-.412 

(.777)- 

(.179)* 

(.363) 

(.403) 

(267) 

km  1 

500 

29.0 

-1.427- 

-.150* 

-.865 

-.918 

-.696 

(2060)- 

(.323)* 

(.855) 

(.938) 

(.644) 

600 

36.9 

-2079- 

-.022* 

-1134 

-1.281 

-1.038 

(4.345)- 

(.571)* 

(1.679) 

(1.800) 

(1.301) 

400 

8.6 

-.070 

.129- 

-.047 

-.053 

-.027* 

(.014)* 

(.056)- 

(.014)* 

(.014)* 

(.014)* 

km  4 

500 

29.1 

-.330 

1.033* 

-.170 

.181 

-.135* 

(.118) 

(1.459)- 

(.051) 

(.055) 

(.043)* 

600 

54.5 

-.853 

4.430* 

-.391 

-.392 

-.356* 

(.734) 

(23.159)* 

(.199) 

(.199) 

(.177)* 

Bathtub 

400 

18.6 

-1.069 

-.185 

-.170 

1.125* 

.063* 

(1.175) 

(.234)* 

(160) 

(1.745)’ 

(.361) 

Pm.  1 

500 

26.1 

-1.722- 

-159 

-102 

1.523 

.046* 

(2.996) 

(.427)* 

(.560) 

(3130)- 

(.608) 

600 

316 

-1452* 

-.362 

-.310 

1.761 

.047* 

(6.043)- 

(.727)* 

(.982) 

(4.490) 

(1.254) 

400 

8.1 

-1.786- 

-1.543 

-1.547 

-.936 

.343* 

(3.218)- 

(2463) 

(1476) 

(1.081) 

(.544)* 

p«  .4 

500 

13.3 

-1370- 

-1.826 

-1.814 

-.825 

.585* 

(5.649)- 

(3.472) 

(3.446) 

(1.031)* 

(1.303) 

600 

18.7 

- 3.072- 

-1191 

-1175 

-.875 

.841* 

(9.466)- 

(5.013) 

(4.983) 

( 1 .285)* 

(2.792) 

*  Bat  estimation  method. 

*  Wont  estimation  method. 
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3.  A  Comparison  of  the  Various  Methods 

A  simulation  siudy  of  data  such  as  that  collected  at  NCTR  was  performed.  Three  groups 
of  48  lifetimes  were  simulated  with  all  testing  stopping  at  280,  420,  and  560  days, 
respectively,  for  the  three  groups.  Distributions  with  mean  survival  times  of  400,  500,  and 
600  days  were  used.  The  generated  lifetimes  greater  than  or  equal  to  the  sacrifice  time  for 
each  particular  group  were  considered  as  censored.  The  remaining  set  of  observed  lifetimes, 
along  with  the  number  censored  at  the  three  sacrifice  times,  constituted  a  single  sample. 
For  each  of  the  distributions  studied,  1000  such  samples  were  generated.  Weibull  distribu¬ 
tions  with  shape  parameters  .5,  decreasing  failure  rate,  1,  constant  failure  rate,  and  4, 

Tibia  2 

Bias/ 100  {and  MSE/tfXf)  for  estimating  90  th  percentile  for  various  methods  of  completion 

Restricted 

Estimated  Weibull  Weibull 
BHK  order  statistic  WTAIL  RWTML 
Distribution  u.  K.-M  extension  extension  extension  extension 
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increasing  failure  rate,  were  used.  Lognormal  distributions,  failure  rate  changes  from 
increasing  to  decreasing,  with  first  two  moments  comparable  to  the  above  Weibull  distri¬ 
butions  with  k  -  1  and  k  ■  4,  were  also  used.  Finally,  a  bathtub  hazard  model  of  Glaser 
(1980),  failure  rate  changes  from  decreasing  to  increasing,  was  used.  This  distribution  is  a 
mixture  of  an  exponential  of  parameter  X  with  probability  1  -  p  and  a  gamma  with 
parameter  X  and  index  3  with  probability  p.  Mixing  parameters  of  p  -  .  1  and  p  ■  .4  were 
used. 

The  bias  and  MSE  for  the  estimation  of  the  tail  probabilities,  i.e.,  the  completed  portion 
of  the  product-limit  estimator,  were  calculated  for  each  hypothesized  distribution  and  for 
each  competing  method  of  completion.  Since  these  results  were  extremely  similar  to  those 
found  in  estimating  mean  survival  time,  m  »  /o  P(t)  dt,  we  show  only  the  bias  and  MSE 
of  each  competing  estimator  of  m  in  Table  1.  This  also  allows  us  to  demonstrate  the 
magnitude  of  the  bias  and  MSE  of  the  product-limit  estimator  of  a.  The  bias  and  MSE  for 
estimating  the  90th  percentile  are  also  presented  for  the  various  estimation  methods  in 
Table  2.  As  one  would  expect,  the  Kaplan-Meier  (K-M)  estimator  performs  considerably 
more  poorly  than  the  other  estimation  schemes.  The  BHK  extension  does  very  well  if  the 
underlying  distribution  is  exponential  or  lognormal  with  Erst  two  moments  compatible 
with  the  exponential.  BHK  does  reasonably  well  for  the  bathtub-shaped  hazard  model,  but 
it  performs  very  poorly  for  the  Weibull  with  increasing  failure  rate  and  for  the  lognormal 
with  first  two  moments  compatible  with  the  WeibulL 

The  remaining  three  extensions  (EOS,  WTAIL,  and  RWTAIL)  appear  to  be  somewhat 
comparable.  Each  of  them  is  best  under  certain  circumstances  although  many  times  the 
biases  and  MSEs  are  so  close  to  one  another  that  they  are  essentially  equivalent  Only  the 
EOS  extension  has  the  desirable  property  of  never  being  worst  It  usually  is  competitive 
with  the  method  that  is  best  Ordering  the  extensions  from  the  standpoint  of  simplicity, 
from  simplest  to  most  complex,  we  have  BHK,  WTAIL,  RWTAIL,  and  EOS. 

In  summary,  the  Kaplan-Meier  estimator  should  probably  be  extended  in  the  presence 
of  extreme  right  censoring.  The  choice  of  extension  depends  on  one’s  knowledge  of  the 
distribution  of  lifetimes  under  consideration  and  the  extent  of  computer  facilities  available. 
If  the  data  follow  an  exponential-type  distribution  or  if  no  computer  facilities  are  present, 
the  BHK  method  is  the  extension  of  choice  due  to  its  simplicity.  If  the  data  exhibit  a 
nonconstant  failure  rate  and  computer  facilities  are  available,  then  the  RWTAIL  or  EOS 
extensions  seem  to  be  advisable. 
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Resume 

On  sait  que  l’estimateur  de  Kaplan-Meier  est  an  estunateur  biaise  de  la  function  de  survie  quand  le 
pourcentage  d’observauons  censurees  est  tres  eleve.  Pluaeurs  modifications  de  I'eshmateur  de  Kaplan- 
Meier  sont  examinees  et  coraparees  du  point  de  vue  de  leun  biais  et  scans  moyens  quadrahques. 
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SUMMARY 

Improved  bounds  on  Che  marginal  survival  function  based  on  data  from  a 
competing  risk  experiment  are  obcained.  These  bounds  are  obtained  by  specifyin 
a  range  of  possible  concordances  for  Che  risks.  These  bounds  are  cighcer  than 
those  of  Peterson  (1976).  A  comparison  Co  ocher  existing  bounds  is  also  made. 


Key  Words:  Competing  Risks, 


Produce  Limic  Estimator.  Ncc  Survival  Funccion, 


Coefficient  of  Concordance. 


I.  INTRODUCTION 


A  common  problem  in  survival  analysis  is  to  estimate  the  marginal 
survival  function  of  the  time,  x,  until  some  event  such  as  remission, 
component  failure,  or  death  due  to  a  specific  cause  occurs.  Often  obser¬ 
vation  of  this  main  event  of  interest  is  impossible  due  to  the  occurrence 
of  a  competing  risk  at  some  time  Y  <  X,  such  as  censoring,  failure  of  a 
different  component  in  a  series  system,  or  death  from  some  cause  not 
related  to  the  study.  Standard  statistical  methods,  which  assume  these 
competing  risks  are  independent,  estimate  the  marginal  survival  function 
by  the  Product  Limit  Estimator  of  Kaplan  and  Meier  (1958)  .  This  estimator 
has  been  shown  to  be  consistent  for  the  marginal  survival  function  by 
Langberg,  Proschan  and  Quinzi  (1981)  when  the  risks  follow  a  constant 
sum  model  defined  by  Williams  and  Lagakos  (1977).  When  the  risks  are 
not  in  the  class  of  constant  sum  models,  the  Product  Limit  Estimator 
is  inconsistent  and,  in  such  cases,  the  investigator  may  be  appreciably 
misled  by  assuming  independence. 
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In  Che  competing  risks  framework  we  observe  T  ■  minimum(x,Y)  and 
I  *  X(X  Y)  where  X(*)  denotes  the  indicator  function.  Tsiatis  (1975) 
and  others  have  shown  that  the  pair  (T,I)  provides  insufficient  information 
co  determine  Che  joint  distribution  of  X  and  Y.  That  is,  there  exists  both 
an  independent  and  a  dependent  model  for  (X,Y)  which  produces  the  same 
joint  distribution  for  (T,I).  However,  these  "equivalent"  independent 
and  dependent  joint  distributions  may  have  quite  different  marginal 
distributions.  Also,  due  to  this  identif iability  problem,  there  may  be 
several  dependent  models  with  different  marginal  structures  which  will 
yield  Che  same  observable  information,  (T,I).  In  light  of  the  consequences 
of  the  uncestable  independence  assumption  in  using  the  Product  Limit 
estimator  to  estimate  the  marginal  survival  funccion  of  X,  it  is  important 
to  have  bounds  on  this  function  based  on  the  observable  random  variables 
(T,I)  and  some  assumptions  on  the  joint  behavior  of  X  and  Y. 

Peterson  (1976)  has  obtained  general  bounds  on  the  marginal,  survival 
funccion  of  X,  S(x),  based  on  the  estimable  joint  distribution  of  (T,I). 

Let  Q^(x)  *  P(T  >  x,  I  ■  1),  and  Q2(x)  *  P(T  >  x,  I  *  0  )  be  the  crude 
survival  functions  of  T.  His  bound,  obtained  from  the  limits  on  the  joint 
distribution  of  (X,  Y)  obtained  by  Frechet  (1951),  is 


Qx(x)  +  Q2(x)  <  S(x)  <  Ql(x)  +  Qx  (0)  .  d-1) 

Since  these  bounds  allow  far  any  dependence  structure,  they  can  be  very  wide 
and  provide  litcle  useful  information  co  an  invest igator. 

Fisher  and  Kanarek  (1974)  have  obtained  tighter  bounds  on  S(x)  in 
terms  of  a  dependence  measure  a.  Their  model  assumes  that  simultaneous 
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co  che  occurrence  of  Y  an  evenc  occurs  which  eicher  scretches  or  concraccs 
Che  remaining  life  of  X  by  an  amounc  associated  wich  a.  ThaC  is, 

P(X  >  xjY  =  y  <  x)  *  P(X  >  y  +  a(x-y)jY  >  y  +  a(x-y)).  A  large  a,  for 
example,  implies  thac  a  small  survival  afcer  censoring  is  che  same  as  a-cimes 
as  much  survival  if  censoring  was  noc  presenc.  They  show  chac  if  a  is  assumed 
known,  chen  che  marginal  survival  function  can  be  estimated  from  the 

a 

observable  information.  Also  these  estimates,  Sa(x),  are  decreasing  in  a. 

For  their  bounds,  the  investigator  specifies  a  range  of  possible  values 
aL  <  a  <  au  so  thac  Sa  (*)  i  S(x)  <_  S  (x). 

U  L 

Recently,  Slud  and  Rubenstein  (1983),  have  proposed  general  bounds. 

They  show  that  knowledge  of  the  function 


P  (x)  -  lira 
6-K) 


P(x  <  X  <  x+6 

X  >  x.  Y  <  x) 

P(x  <  X  <  x  +<5 

X  >  x,  Y  >  x) 

along  with  the  observable  information  (T,I)  is  sufficient  to  uniquely 
determine  the  marginal  distribution  of  X.  These  estimates  Sp(t)  are 

decreasing  functions  of  p  for  fixed  x.  Their  bounds  are  obtained  by 
specifying  a  range  of  possible  values  p^(x)  <_  p(x)  <_  P2<x)  so  that  if 
o(x)  is  che  true  function  S^Cx)  £  S(x)  <_  Sp^(x). 

In  this  paper  we  obtain  different  bounds  on  the  marginal  survival 
function  by  assuming  a  particular  dependence  scructure  on  X  and  Y.  These 
bounds  are  functions  of  the  observables  (T,I)  and  a  familiar  dependence 
measure,  che  concordance  probability  between  X  and  Y.  In  Section  2  we 
describe  this  model  in  detail.  In  Section  3  we  derive  Che  bounds  and  show 
consiscency  when  che  dependence  parameter  is  known.  In  seccion  U  these 
bounds  are  compared  to  those  obtained  by  Peterson,  Fisher  and  Kanarek, 


and  Slud  and  Rubenstein. 


u 

II.  THE  MODEL 

The  dependence  structure  we  shall  employ  to  model  the  joint  survival 
was  first  introduced  by  Clayton  (1978)  to  model  association  in  bivariate 
LifetableS  and,  later,  by  Oakes  (1982)  to  model  bivariate  survival  data. 
Let  S(x)  *  P(X  >  x),  R(y)  -  P(Y  >  y),with  S(0)  -  R(0)  -  1,  be  the 
continuous  univariate  survival  functions  of  the  death  and  censoring  times, 
respectively.  For  9  _>  1  define  F(x,y)  -  P(X  >  x,  Y  >  y)  by 


F(x,y) 


({ 


1  j9"1 


S(x) 


r  1  ;  9-1  ,  ,-1/ (9—1) 

W}  _1  1 


(2.1) 


This  joint  distribution  has  marginals  S  and  R.  As  9-*-l,  then  (2.1)  reduces 
to  the  joint  distribution  with  independent  marginals.  For  F(x,y)  -*• .. 
min(S(x),  R (y ) )  the  bivariate  distribution  with  maximal  positive  association 
for  these  marginals.  The  probability  of  concordance  is  9/(9  +  1)  so  that 
Kendall's  (1962)  coefficient  of  concordance  is  T  *  (8  -  l)/(9  +  1)  which 
spans  Che  range  0  to  1. 

This  model  has  a  nice  physical  interprecation  in  terms  of  the 
functions  A(x|Y  *  y)  and  A(x|Y  >  y) ,  the  hazard  functions  of  X  given  Y  «  y 
and  X  given  Y  >  y,  respectively.  From  (2.1)  one  can  show  that 

l(x|Y  -  y)  3  9A(x|Y  >  y) 
or 

P(X  >  x|Y  -  y)  -  (P(X  >  x | Y  >  y ) ] ?  (2.2) 

For  9  >  1  the  hazard  rate  of  survival  if  censoring  occurs  at  time  y  is 
9  times  the  hazard  rate  of  survival  if  censoring  does  not  occur  at 


time  y.  This  implies  that  the  hazard  rate  after  censoring  occurs  is 
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accelerated  by  a  factor  of  6  over  che  hazard  rate  if  censoring  had  noc 
occurred.  Also  when  9  ■  1 ,  (2.2),  reduces  to  the  condition  required  by 
Williams  and  Lagakos  (1977)  for  a  model  to  be  constant  sum  and  hence  for 
the  usual  product  limit  estimator  of  S(t)  to  be  consistent  (See  Basu  and 
Klein  (1982)  for  details). 

Oakes  (1982)  also  shows  that  (2.1)  can  be  obtained  from  the  following 
random  effects  model.  Let  S*(x)  *  exp  {-  ^s~('x)~^  ^  an(*  ^eC  be 

similarly  defined.  Let  W  have  a  gamma  distribution  with  density 
-1  -w 

g(w)aw°  e  and  conditional  on  W  ■  w  let  X,Y  be  independent  with 

survival  functions  (S*(x)  }w  and  {R*(y)}w.  Then,  unconditionally,  X,Y  have 
the  joint  survival  function  F(x,y)  given  by  (2.1). 

For  fixed  marginals  S  and  R  the  joint  probability  density  function, 
f(x,y),  can  be  shown  to  be  totally  positive  of  order  2  for  all  8  >_  1. 

This  implies  that  (X,Y)  are  positive  quadrant  dependent.  In  particular, 
one  can  show  that  for  S,R  fixed  the  family  of  distributions 
F  ■  (F(x,y):  9^1}  is  increasing  positive  quadrant  dependent  in  9  as 

defined  by  Ahmed,  et  al .  (1979). 
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III.  80UNDS  ON  MARCINAL  SURVIVAL 

Suppose  chac  X  and  Y  have  Che  joinC  distribution  (2.1)  and  let 
T  -  min  (X,Y),  then  the  survival  function  of  T  is 


,  .  , ,  1  ,9-1  ,  1  , 9-1  , , -  t — r 

F(T)  ’  ^S(t)  *  R(t)  "  ^  9-1 


(3.1) 


and  the  crude  density  function  associated  with  X, 


q,(c)  *  77  p(T  <  c,  X  <  Y),  is  given  by 


q,(c)  -  . 

S9(t) 


(3.2) 


where  s(t)  *  -dS(t)/dt. 

Now  consider  the  differential  equation 


s(t)/S9(t)  -  qi(tyiF(c)f 


(3.3) 


and  suppose  0  is  known.  Then  the  solution  of  (3.3)  for  S(t)  is 


’  1  *  ^  i  du] 


(3.4) 


t  q,(u) 

exp<  '  '  75)  d'J> 


if  0  -  1 , 


The  functions  F(*)  and  q^(*)  are  directly  estimable  from  the  data  one 
sees  in  a  competing  risks  experiment.  Let  T^ ,  ....  Tr  denote  the  observed 
test  times  of  n  individuals  put  on  test  and  let  1^,  i  *  1,  ....  n  be  1  or  0 
according  to  whether  the  T^  was  an  observation  on  or  Y^,  respectively. 


jOw  SU 
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n  n 

Define  K(t )  *  lx  (T .  >  c)/n  and  Q,  ( t )  *  E  x  (T .  £  C ,  L.  *  1). 

i-1  1  1  1-1  1  1 


Then  if  9  is  known. 


a  natural  estimator  of  S  (t)  is 

u 


c  l 

S0(t)  -  [1+  (0-1)  £  [?(u)p  1"  (0-1) 


t  dQ. (u) 

-  exp(  -  I  -  ) 

0  F(u) 


if  9  >  1 


if  0  -  1 


(3.5) 


For  0-1,  this  estimator  is  of  the  form  of  the  hazard  rate  estimator 
proposed  by  Nelson  (1972).  The  estimators  (3.5)  can  be  expressed  in 
the  following  form  for  computation  purposes. 


SQ(t) 


exp 


(n-i+1) 


T  <t  I 

(i)-  i’  (i) 


if  9 


[1  +  (8-l)n 


9-1 


I  (n-i+1) 


"  (PI 


T(i)<t,I(i)"1 


(3.6) 


if  9  >  1 


where  T,,,,  ....  T,  .  are  the  ordered  death  times, 
in; 

For  9  known  and  if  the  true  underlying  joint  distribution  of  (X,Y) 


is  of  the  form  (2.1)  then  S0(t)  is  a  consistent  estimator  of  S(t)  as  shown 
by  the  following  theorem. 


Theorem  1.  Let  (X,Y)  have  the  form  (2.1)  with  marginals  S(t),  R(t) 
respectively.  Let  9  M  be  known.  Then  on  the  set  where  S(t)  >  0  we  have 
SQ(t)  -  S(t)  a.s. 


6 

Proof : 

For  9  »  1,  che  result  follows  by  a  theorem  of  Langberg,  Proschan  and 
Quimi  (1981).  Suppose  that  9  >  1.  Note  that  Q^(t)  ■*  Q^(t)  a.s.  and 

A  A 

F  (u)  ■*  F(u)  a.s.  by  the  strong  law  of  large  numbers.  Since  SQ(c)  is  a 

r  0 

dQl (u) 

continuous  function  of  *  — — g  in  the  support  of  F(u),  it  suffices  to  show 

0  [ F(u) J 

dQ, (u)  C  dQ, (u) 

I  -1 -  -  t  —I _  a.s. 

0  (F(u)]0  0  [F(u)j9 

Now,  after  an  integration  by  parts. 


dQ  (u)  Q  (C)  „ 

‘  - - Q  *  7 - q-  '  ) 

0  [F(u)]°  (F(t ) ]  0  F°(u) 


QiU)  C  -  ,  C  . 

1  - - 5-  ,  (Q.(u)  -  Q  (u)Jd(i  )  +  j  Q.  (u)d(i  ) 

[F(t)]9  Q  1  1  F9(u)  q  1  F0(u) 


Q  (t)  -  Q  .(t)  „ 

'  7 - z - '  (Q,(u)  -  Q.  (u)]d(i  ) 

[F  (u)]a  0  1  1  F°(u) 


C  dQ, (u) 

♦  '  ^-9  • 

0  F  (u)° 

By  the  dominated  convergence  theorem 


(3.7) 


i  im 
n-*» 


C  dQ, (u) 
f  1 

0  [F(u)f 


t 

( 

0 


dQj^  (u) 
(F(u)f 


a.s., 
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11.  V‘>  -  V‘> 


n-*® 


[F  (u)]9 


0.'  a.s.. 


and 


lira  sup  {jQ^Cu)  -  (u) | }  -  0,  a.s. 


n-*® 


Hence,  applying  the  above  results  to  (3.7),  the  result  now  follows:  // 

To  obtain  bounds  on  Che  net  survival  function  based  on  data  from  a 

competing  risks  experiment,  we  proceed  as  follows.  First,  note  chat  from 

(3.5)  it  is  true  that  §  (t)  is  a  decreasing  function  of  0  for  fixed  t. 

£ 

Also,  as  0  -*•  1+  we  have  Sg(t)  t  exp  (  —  /  F  *(u)dQ^ (u)) . 


which  provides  an  upper  bound.  Notice  that  this  upper  bound  corresponds 

A  A 

to  an  assumption  of  independence.  As  9  ■*  ®  one  can  show  that  Sg(t)  1  F(t) 

which  corresponds  to  Peterson's  (1976)  lower  bound. 

# 

la  practice  the  above  bounds,  with  9*1,®,  while  shorter  than 
Peterson's  bounds,  may  still  be  quite  wide. 

Tighter  bounds  may  be  obtained  by  an  investigator  specifying 
a  range  of  possible  values  for  9.  If  the  sample  size  is  sufficiently 
large  and  <  0  <_  02,  Chen  S0  (t)  <  S(t)  <_  Sg  (t).  Specifying  0^  92 

is  equivalent  to  specifying  a  range  of  values  <  x  <  for  the 
coefficient  of  concordance  t  since  9  ■  (1 +t)/(1  --T) .  Hence  the  primary 
value  of  S 0 ( c )  is  in  putting  bounds  on  S(t)  rather  chan  on  estimation  of 
S(t). 
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IV.  EXAMPLE  AND  COMPARISONS 

To  illustrate  the  bounds  obtained  in  the  previous  section,  consider 
the  mortality  data  reported  in  Hoel  (1972).  The  data  was  collected  on  a 
group  of  RFM  strain  male  mice  who  were  subjected  to  a  dose  of  300  rads  of 
radiation  at  age  5-6  weeks.  There  were  three  competing  risks,  thymic 
lymphoma,  reticulutum  cell  sarcoma,  and  ocher  causes  of  death.  For 
illustrative  purposes  we  consider  reticulum  cell  sarcoma  as  Che  risk  of 
interest . 

Table  1  reports  the  value  of  Sfl(t)  for  concordance  r  ■  (8  -  1 ) / (0  +  1) 

a 

The  value  of  Sg(c)  at  T  ■  0  corresponds  to  Nelson's  (1972)  hazard  rate 
estimator  assuming  independence.  Peterson's  upper  and  lower  bound 
(t  *  1)  are  also  reported  as  are  Fisher  and  Kanerek's  bounds  and  the  Slud 
and  Rubens tein  bounds  for  several  values  which  reflect  a  positive 
association  between  risks. 

From  Table  1  we  first  note  that  Peterson's  bounds  are  very  wide. 
Substantial  improvement  is  obtained  if  one  assumes  a  non-negative 
dependence  structure  between  risks  (See  Table  2).  Further  tightening 
of  these  bounds  is  achieved  by  assuming  chat  T  is  in  the  range  0  to  .5 
where  the  width  of  the  boundaries  is  at  most  about  50Z  of  that  of  Peterson' 
bounds. 

Substantial  improvement  in  the  general  bounds  is  also  obtained  by 
the  bounds  of  Fisher  and  Kanerek  or  Slud  and  Rubenstein.  The  bounds  of 
Fisher  and  Kanerek  assume  a  specific  censoring  pattern  and  require  a 
specification  of  a  stretching  constant  a.  Without  some  additional  informa¬ 
tion,  such  specification  may  be  impossible.  Slud  and  Rubenstein's  bounds 
are  for  the  general  dependence  structure.  Their  bounds  require  the 
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specification  of  Che  o(t)  function.  This  function  is  a  quantity  which 
is  not  easily  conceptualized  by  investigators  from  either  a  scacistical 
or  biological  perspective.  This  makes  it  questionable  whether  reasoable 
upper  and  lower  bounds  on  p(t)  can  be  extracted  from  one's  prior  beliefs. 
The  major  advantage  of  the  bounds  printed  in  this  paper  is  that  they 
require  only  the  specification  of  an  upper  and  lower  concordance,  a 
measure  quite  familiar  to  most  investigators  and  easily  explainable  to 
nonscat is ticians. 


t 


( 

a 

r> 
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Table  2 

RELATIVE  SIZE  OF  THE  BOUNDS  ON  NET  SURVIVAL 
FOR  AN  ASSUMED  DEPENDENCE  STRUCTURE 
AS  COMPARED  TO  PETERSON'S  BOUNDS 


Time 

0  <  T  <  1 

0  <_  T  <_  .5 

0  1  T 

350 

.9707 

.0879 

.2674 

525 

.9352 

.2449 

.5931 

600 

.7338 

.5171 

.6787 

620 

.6722 

.5120 

.6298 

650 

.5009 

.4420 

.4870 

675 

.3831 

.3576 

.3797 

700 

.2883 

.2767 

.2833 

750 

.0600 

.0600 

.0600 
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1.  Introduction 


Consider  a  system  consisting  of  several  components  linked  in 
series.  For  such  a  system  the  failure  of  any  one  of  the  components 
causes  the  system  ho  fail.  Ln  a  biological  or  medical  context 
we  can  consider  the  components  to  be  different  lethal  diseases 
and/or  different  reasons  for  removal  from  the  study.  In  a  clini¬ 
cal  trials  framework  the  primary  response  of  interest,  death  or 
remission,  and  censoring  can  be  considered  as  components  of  the 
system.  This  general  formulation  has  been  detailed  in  the  theory 
of  competing  risks  (cf.  David  and  Moeschberger  (1978)). 

A  cannon  assumption  in  such  a  formulation  is  that  the 
component  lifetimes  are  statistically  independent.  Several 
authors  have  shown  that  based  on  data  frcm  series  systems  only , 
this  assumption,  by  itself,  is  not  testable  because  there  is  no 
way  to  distinguish  between  independent  or  dependent  component 
lifetimes  (see  Basu  (1981),  Basu  and  Klein  (1982),  Miller  (1977), 
Peterson  (1976),  etc.).  However,  several  authors  (see  Lagakos 
(1979)  p.  152  and  Easterling  (1980)  p.  131)  have  pointed  out  the 
need  to  determine,  quantitatively,  how  far  off  one  might  be  if  an 
analysis  is  based  on  an  incorrect  assumption  of  independence. 

To  study  the  effects  of  erroneously  assuming  independence  we 
shall  assume  that  each  of  the  components  have  exponentially  dis¬ 
tributed  lifetimes  when  tested  separately  and  that  the  property  of 
marginal  exponent iality  will  be  preserved  even  though  some 
dependence  nay  be  induced  when  the  components  are  linked  in  series. 
The  assumption  of  exponentially  distributed  component  lifetimes 
has  been  made  by  Mann  and  Grubbs  (1974) ,  when  finding  confidence 
bounds  on  system  reliability,  Boardman  and  Kendall  (1970),  Uien 
estimating  component  lifetimes  frcm  system  data,  and  Miyaraura 
(1982),  vften  combining  component  and  sytem  data.  (See  Barlow  and 
Proschan  0.975)  or  Mann,  Schaffer,  and  Singpurwalla  (1974)  for  a 


acre  complete  review.}  We  shall  model  the  dependence  structure  by 
the  three  models  of  Gumhel  (I960),  a  model  proposed  by  Downton 
C197Q) ,  and  a  model  described  by  Oakes  (1982).  These  models  are 
briefly  described  in  Section  2. 

The  effects  of  a  departure  from  the  assumption  of  independent 
component  lifetimes  will  be  addressed  for  two  distinct  situations. 
The  first  situation  arises  in  modeling  the  performance  of  a 
theoretical  series  system  constructed  from  two  components.  Here, 
based  on  testing  each  component  separately  or  on  engineering 
design  principles,  it  is  reasonable  to  assume  that  the  components 
are  exponentially  distributed  with  known  parameter  values .  Based 
on  this  information,  we  wish  to  predict  parameters  such  as  the 
mean  life  or  reliability  of  a  series  system  constructed  from 
these  components.  In  Section  3  we  describe  how  these  quantities 
are  affected  by  departures  from  independence. 

The  second  situation  involves  making  inferences  about 
component  lifetime  distributions  from  data  collected  on  series 
systems.  Conmonly,  data  collected  on  such  systems  are  analyzed 
by  assuming  a  constant-sun  model,  of  which  independence  is  a 
special  case  (compare  Williams  and  Lagakos  (1977)  and  Lagakos  and 
Williams  (1978)).  In  Section  4.1  we  study  the  properties  of  the 
maximum  likelihood  estimators  of  the  component  mean  life  calcu¬ 
lated  under  an  erroneous  assumption  of  independent  exponential 
component  lifetimes  as  mentioned  above.  Because  of  the  wide 
spread  use  of  the  ncn parametric  estimator  proposed  by  Kaplan- 
Meier  (1958)  for  the  component  reliability  we  study  in  Section  4.2 
its  properties,  when  the  marginal  reliahil  ities  are  exponential 
and  independence  is  incorrectly  assumed. 

2.  The  Models 

Consider  a  two  component  series  system  with  component  life 
lengths  X^,  X^.  Suppose  that  each  X^  has  an  exponential  survival 
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function 

(2.1)  f\.Ct)  a  P(X£  >  t)  =  expC-^  t),  Ai  >  0,  t  >  Q. 

This  assumption  is  trade  on  the  basis  of  extensive  testing  of 
each  component  separately  or  on  knowledge  of  the  underlying 
mechanism  of  failure. 

To  examine  the  effects  of  a  departure  from  independence  we 
consider  five  bivariate  exponential  models,  each  with  marginals 
equivalent  to  (2.1) .  The  first  three  models  are  due  to  Gumbel 
(1960);  the  last  two  models  are  due  to  Down  ton  (1970)  and  Oakes 
(1982). 

2.1  Gumbel 's  Model  A 

For  this  model  the  joint  survival  function  is 

(2.2)  P(X1  >  Xj  >  Xj)  =  exp(-  X^  -  X^  -  ^12xix2)» 

x]>»  *2  —  °»  *1»  *2  >  °»  0  i  *12  -  *1*2* 

The  correlation  between  X.  ,  X-  is 
XX  A  c 

p  =  -  -i-i  exp(X1X2/X12)  E^-  ^1^2/*12)  " 

12  <■ 

where  E- (z)  =  f*  exp(-u)  du  is  the  integrated  logarithm 

u 

For  this  model  o  varies  from  -  .40365  to  0  as  X^2  decreases  from 
X^X2  to  0.  It  is  never  positive.  The  regression  X2  on  X^  is  non- 
leanear  with 

2 

E(X-  jX2  =  x2)  =  (XL  ♦  X12x2  -  X12/X2)/(X1  ♦  X12x2)  . 

2.2  Gunbel’s  Model  B 

For  this  model  the  joint  survival  function  is 

(2.3)  P(Xj>  x^  >  x2)  =  exp(-  X^  -  X2x2){l  ♦  4p(l-exp(-X1x1) ) 

Q-  exp(-  X^)},  X1#  X2  >'0»  ^  Xj  v  0,  -  1/4  <  p  <  1/4. 
The  correlation,  P,  nay  be  positive  or  negative.  The  regression 


of  X2  on  X^  is  again  nonlinear  with. 

EO^|X2  a  x2)  a  {l  f  2P  -  >40  expC-  X2x2))A1- 

The  effects  of  a  departure  from  independence  on  modeling  system 
reliability  and  estimating  component  reliabilities  has  been  studied 
in  detail  in  Moeschberger  and  Klein  (1984). 

2.3  Gumbel’s  Model  C 


For  this  model  the  joint  survival  function  is 


P(XX  >  V  *2  >  V  =  exP{“  CCX1xl)ni  *  (X2x2)m]1/m}  ’ 


^1’  X2  >  °’  m  -  *1’  *2  >  °' 


The  correlation  is 


P  =  (4 


*  2m)  ‘ 


(cos  0  sin  9  )“ 

I  ~  ■  2+2/m  de  -  l 
0  faos  8®  ♦  sin  0°) 


which  varies  frcm  0  to  1.  For  this  model  m  =  1  corresponds  to 
independence  and  as  m  ♦  - 

(2.S)  P(X^  >  5Cj_,  ^  >  *2^  *  minimum  [exp(-  IjX^),  exp(-  1^)] 
the  Frechet  (1958)  upper  bound  for  these  marginals. 

2.4  Downton's  Model 

Downton  (1970)  suggests  modeling  bivariate  exponential  systems 
by  a  successive  damage  model .  This  model  assumes  that  in  a  tvo 
component  system  the  times  between  successive  shocks  on  each 
component  have  independent  exponential  distributions-  and  that  the 
nunber  of  shocks  required  to  cause  each  component  to  fail  follows 
a  bivariate  geometric  distribution.  The  joint  probability 
density  function  of  the  component  life  times  is 


(2.5)  f(x1',x2)  =  exp/-  (X1x1  .♦  X2x2) 


where  IqC*)  is  the  modified  Bessel  function  of  the  first  kind  of 
order  zero,  and  X^,  X2  >  Q»  x^,  Q,  Q  <  p  ^  1.  The  correla¬ 

tion  between  X^,  ^  is  p  which  spans  the  interval  [Q,l].  As  P  *  1 
the  joint  survival  function  of  apprrvkohes  the  upper 

Frechet  distribution  (2.S).  For  this  model 
ecx1|x2  =  5^}  =  ci-p/Xj^  *p  x2X2/xr 

4 

2.S  Oakes'  Model 

Oakes  (1982)  has  proposed  a  model  for  bivariate  survival 
data.  This  model  was  first  proposed  by  Clayton  (1978)  to  model 
association  in  bivariate  lifetables.  Special  cases  of  Oakes’ 
general  model  have  been  suggested  by  Lindley and Singpurualla  (1985) 
and  Hutchinson  0.981). 

For  this  model  the  joint  survival  probability  is 

(2.7)  P(X1>  x1,  X2>  x? )  =  Cexp(X1(8-l)x1)  ♦  «p(X2<e-l>x2>-l3“1/d~*> 

where  X^,  X 2  >  0,  8^1,  x^,  X2  ^  0. 

For  9  si,  Xj,,  X2  are  independent  and  PO^  >  xi»  x2  >  x^  ♦  (2.5) 
as  8  For  this  model  Kendall's  (1962)  coeffcient  of  concordance 

is  t  =  (9-l)/(8+l)  which  spans  the  range  0  to  1.  The  correlation, 

P,  also  spans  the  range  0  to  1  and  is  found  numerically. 

This  model  has  the  following  physical  interpretation.  Let 
r(XjJX2  =  x^)  and  r(x^|X2  >  x^)  be  the  conditional  failure  rates 
of  X^  given  =  x^  and  X2  >  x^  Then  Kx^|X2  -  x^  -  9  r(x^|X2  >  x^) 

The  model  can  also  be  derived  from  a  random  effects  model. 

This  formulation  assumes  that  when  the  conponents  are  tested 
separately  under  ideal  conditions  the  component  survival  functions 
are  S^(t)  =  expC-  exp(X^t(9-l) )  *  1],  i  =  1,  2,  and  that  when  the 
two  components  are  put  in  a  series  system  in  the  operating 
environment  there  is  a  random  factor  W  which  simultaneously  changes 
each  ccnponent  life  distribution  to  S^(t) .,  jL£  W  has  a  gpnina 

distribution  with  density  function  f  (x)  aw^"^  e-**  then, 

-S- 


unconditionally,  the  joint  survival  function  (2.7)  holds. 
2.6  Fr4cfiet  Sounds 


Frechet  (1958)  obtained  bounds  of  the  joint  survival  functions 
fchich  can  be  obtained  for  any  set  of  marginal  distributions.  For 
exponential  marginals  these  are 


-\.x. 

MAXIMUM  (e  1  *  e  -  1,  0)  £  P^  >  5^, 


-Xlxl 


(2  >  X2)  i 

Q^2*2) 


MINIMUM  (e 

For  this  set  of  marginals  the  lower  Frechet  distribution  has 
correlation  -.694  and  the  upper  Frechet  distribution  has 
correlation  1.0.  These  are  the  minimal  and  maximal  correlations 
for  exponential  marginals. 


3.  Errors  in  Modeling  System  Life 

Suppose  that  based  on  extensive  testing  or  based  on  theore¬ 
tical  considerations  each  of  the  two  components  in  a  series 
system  is  known  to  have  an  exponential  distribution,  (2.1)  with 
marginal  means  1/X^,  l/X^,  respectively.  It  is  of  interest  to 
predict  the  system  reliability  F(t)  =  P(X^  >  t,  X^  >  t)  and  the 

system  mean  life  u  =  FCt)  dt.  If  the  investigator  assumes  that 
the  two  components  are  independent  then  the  system  reliability  is 

(3.1)  Fj(t)  =  exp(  -U1  +  X2)t>  and  system  mean  life  is  Uj  =  l/(*^  *  x 

If  the  components  are  not  independent,  but  in  fact  follow  one 
of  the  models  in  Section  2,  then  a  measure  of  the  effects  of 
incorrectly  assuming  independence  'A(t)  -  (T(t)  -  Fj(t))/Fj(t)  and 
d  =  (u  -  Wl)/ur  for  predicting  system  reliability  and  system  mean 
life,  respectively,  vhere  fit)  and  u  are  computed  under  the 
appropriate  dependent  model.  Values  of  F(t)  can  be  computed 
directly  fran  (2.2),  (2.3),  (2.4),  (2.6)  (by  nunerical  integration) 
or  from  (2.7).  Expressions  f or  u  are  given  in  Appendix  1.  All 


expressions  for  A(t )  and  $  depend  on  the  values  of  X^  and  X2  only 
thrxxigh  the  ratio  X^/X2  a  K  and  for  K  <  1  the  values  are  equiva¬ 
lent  to  those  for  r  :  For  the  upper  Prechet  distribution, 

K 

-  K5T 

ACtp)  a  p  -  1  vhere  K  £  1  and  tp  =  the  point 

vhere  F^-Ct^)  =  p.  Also  5  =  1/K  for  K  >  1.  For  the  lower  Frechet 

distribution 


-IK  K  1 

A(tp)  =  p“  *  p“  K+l  -  p-1  -  1  if  pK  ^  ♦  _  l  >  0 


otherwise 


and  5  =  K2  +  K  »  1  -  GC+1)?  Y  ♦  (K+l)  4n(Y)  vhere  Y  is  the  solu¬ 

te  K 

tion  of  the  equality  /  ♦  X  s  1.  Table  1  gives  the  values  of 
A(tp)  x  100%  and  <5  x  100%  for  p  s  .9,  .7,  .5,  .3,  -1  for  the  upper 
and  lower  Fr&het  distributions. 

From  Table  1  we  see  that  the  Largest  percent  error  occurs 
when  the  parameters  are  equal.  Also  for  fixed  K  there  is  rela¬ 
tively  snail  error  in  estimating  system  reliability  by  modeling 
a  dependent  system  by  an  independent  system  then  F(t)  is  large. 

For  smaller  values  of  system  reliability  one  can  be  appreciably 
misled.  Errors  in  estimating  system  mean  life  appear  to  be 
substantial  unless  one  component  has  considerably  longer  marginal 
life  than  the  second  one.  In  that  instance,  one  can  see  instinc¬ 
tively  that  the  correlation  would  have  a  minimal  impact. 

Figures  1A-1C  and  2A-2C  are  plots  of  A(tp)  for  p  s  .2Hr  .5,  .75, 
X^  =  1  and  X2  =  1,  1.5  for  the  five  models  described  in  Section  2. 
Figures  3A,  3B  are  plots  of  <S  for,  all  five  models  as  a.  function  of 
the  correlation.  From  these  plots-  note  that  for  positive  correla¬ 
tions  the  most  modeling  error  occurs  for  the  GXmfcel  C  model.  For 
relatively  small  correlation,  .-.2S-  <,  p  <  .25  there  may  still  be  a 
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moderate  Modeling  error,  on  the  order  at  least  ^  1Q%  for  pre¬ 
dicting  system  reliability  at 
estimating  the  mean  system  life. 


FjCt)  =  .2S,  or  FjCt)  =  .S  and  for 


TABLE  t 

UPPER  AND  LOWER  MUX0S  OR  THE  PERCENT  ERROR  IN  WDELIN6  STSTER  LIFE 


FUI-0.9 


F(T)«0.7 


'  Ftri-o.s 


F1T)»0.2 


F<T>«0.1 


HE AN  LIFE 


t 

LONER** 

80UN0 

’upper 

BOUND 

LOWER 

B0UN0 

UPPER 

BOUND 

loner" 

BOUND 

VkR 

BOUND 

LONER 

BOUND 

UPPER  LONER 
BOUND  BOUND 

UPPER 

BOUND 

LONER 

8QUN0 

UPPER 

BOUND 

T~ 

-0.29 

3.41 

-3.81 

19.32 

-17.1* 

4t.  42 

-100.00 

82.37  -100.00  214.23 

-38.43  100.00 

2 

-0.2* 

3.37 

-3.39 

12.42 

-13.27 

23.99 

-100.00 

49.38  -100.00 

113.44 

-37.07 

30.00 

3 

-0.22 

2.47 

-2.  a* 

9.33 

-12.90 

18.92 

-31.32 

33.12  -100.00 

77.83 

-34.83 

33.33 

4 

-0.19 

2.13 

-2.44 

7.39 

-11.02 

14.87 

-44.11 

27.23  -100.00 

38.49 

-32.83 

23.00 

3 

-0.14 

1.77 

-2.12 

4.12 

-9.37 

12.23 

-38.38 

22.22  -100.00 

44.78 

-31.0* 

20.00 

A 

-0.14 

1.32 

-1.87 

3.23 

-8.43 

10.41 

-33.91 

18.77  -100.00 

38.93 

-29.33 

14.47 

7 

-0.13 

1.33 

-1.47 

4.34 

-7.33 

9.03 

-30.33 

14.24  -tOO.OO 

33.33 

-28.18 

14.29 

8 

-0.12 

1.18 

-1.51 

4.04 

-4.82 

8.01 

-27.42 

14.31  -100.00 

29.13 

-26.99 

12.30 

9 

-0.11 

1.04 

-1.37 

3.43 

-4.22 

7.18 

-23.02 

12.79  -100.00 

23.89 

-23.93 

11. 11 

10 

-0.10 

0.94 

-1.24 

3.30 

-3.71 

4.30 

-22.99 

11.37  -100.00 

23.28 

-24.98 

10.00 

11 

-0.09 

o.aa 

-1.17 

3.02 

-5.28 

3.93 

-21.27 

10.33  -100.00 

21.15 

-24.12 

9.09 

12 

-0.08 

0.81 

-1.08 

2.78 

-4.91 

5.48 

-19.78 

9.70  -100.00 

19.38 

-23.34 

8.33 

13 

-o.oa 

0.74 

-1. 01 

2.38 

-4.39 

5.08 

.  -18.49 

8.98  -100.00 

17.88 

-22.42 

7.49 

l< 

-0.07 

0.70 

-0.93 

2.41 

-4.30 

4.73 

-17.33 

8.34  -100.00 

14.59 

-21.94 

7.14 

13 

-0.07 

0.44 

-0.90 

2.23 

-4.03 

4.43 

-14.33 

7.82  -100.00 

13.48 

-21.33 

4.47 

14 

-0.04 

0.42 

-0.83 

2.12 

-3.83 

4.14 

-13.43 

7.34  -100.00 

14.30 

-20.79 

4.23 

17 

-0.04 

0.39 

-0.80 

2.00 

-3.43 

3.93 

-14.43 

4.92  -100.00 

13.43 

-20.24 

3.98 

ia 

-0.04 

0.34 

-0.74 

1.89 

-3.43 

3.72 

-13.93 

4.34  -100.00 

12.88 

-19.74 

5.34 

19 

-0.04 

0.33 

*0.73 

1.80 

-3.29 

3.33 

-13.27 

4.20  -100.00 

12.20 

-19.30 

5.24 

20 

-0.03 

0.30 

-0.49 

1.71 

-3.14 

3.34 

-12.47 

3.90  -92.24 

11.39 

-18.87 

3.00 

21 

-0.03 

0.41 

-0.4* 

1.43 

-3.00 

3.20 

-12.13 

3.43  -88.34 

11.03 

-18.44 

4.74 

22 

-0.03 

0.4* 

•0.44 

1.3* 

-2.88 

3.0* 

-11.43 

3.37  -84. 73 

10.33 

-18.08 

4.33 

23 

-0.03 

0.44 

-0.41 

1.30 

-2.7* 

2.93 

-11.14 

3.14  -81.41 

10.07 

-17.71 

4.33 

24 

-0.04 

0.42 

-0.39 

1.44 

-2.4* 

2.81 

-10.74 

4.93  -78.34 

9.43 

-17.37 

4.17 

23 

-0.04 

0.41 

-0.37 

1.38 

-2.34 

2.70 

-10.34 

4.74  -73.49 

9.24 

-17.04 

4.00 
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4.  Errors  in  Estimating  Component  Parameters 
4.1  Parametric  Estimation 

In  this  section  we  examine  the  effects  of  incorrectly  assuming 
independence  an  the  magnitude  of  the  estimation  error  in  esti¬ 
mating  the  first  component  mean  life  based  on  data  from  series 
systems.  Suppose  that  n  series  systems  are  put  on  test.  For 
each  system  we  observe  the  system  failure  time  and  which  component 
caused  the  failure.  Let  n -  dentoe  the  number  of  systems  where 
the  system  failure  was  caused  by  failure  of  the  im  component , 
i  =  1,  2,  and  let  T  be  the  total  time  on  test  for  all  n  systems. 

If  we  assume  that  the  component  lifetimes  are  independent  and 
exponential  1  y  distributed  then  Moeschberger  and  David  (1971)  show 
that  the  maximum  likelihood  estimator  of  u^,  the  first  component 
mean  life  is 

a 

(4.1)  =  T/h.^  for  >  0. 

This  estimator  is  asymptotically  unbiased  and  for  n  finite 
E(u^)  =  E(T) -Ed/n^i >  0)  due  to  the  independence  of  T  and  n^.. 

Suppose  now  that  the  two  component  lifetimes  are  not  inde¬ 
pendent  but  follow  one  of  the  models  discribed  in  Section  2.  If 
we  incorrectly  assume  independence  then  a  measure  of  the  excess 
bias  due  to  incorrectly  assuming  independence  is 
3  =  [E(u^  |  Dependent  model)  -  E(u^|  independence)  ]/u^.  For  each  of 
the  dependent  models  under  consideration  T  and  are  independent. 
For  large  n,  B  converges  to  (u/p  -  /u-j_  where  u  is  the  mean 

system  life  and  p  is  the  probability  the  first  component  fails 
first,  computed  under  the  dependent  model.  For  finite  n, 

E(u^)  =  n  u  EpO/njJn^  >  Q)  computed  under  the  appropriate  model, 

where  EpC/nJ^  >  Q)  =  £  (£)pKCl-p)K/K  /  (l-(l-p)n).  Expressions 

K=  1 

for  u  and  p  are  given  in  Appendix  1  and  Appendix  2,  respectively. 
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The  expressions  depend  on  A^,  Aj  only  through  the  ratio  K  = 

For  all  models,  p  a  1/2  viien  K  =  1. 

For  the  upper  Fr&het  distribution  p  =  Q  if  K  <  1;  1/2  if  K 
and  1  if  K  >  1.  Hence  for  K  <  1  no  failures  from  the  first  com¬ 
ponent  are  ever  observed  so  that  the  modeling  error  B  becomes 
infinite  for  all  n.  For  K  >  1,  p  =  1  and  p  =  so  that 
B  =  Cl  -  ECu^l  Independence) /u^)  vfiich  tends  to  0  as  n  ♦  In 
this  case  the  models  with  correlation  ranging  from  Q  to  1  have  B 
increasing  for  p  <  pQ  and  decreasing  for  p  >  pQ.  For  the  lower 
Frechet  distribution,  p  is  the  value  of  X  uhich  solves  the 
equation  X  ♦  X1^  -  1  =  0.  For  K  <  1  we  have  p  <  1/2  and  for  K  >  1 
we  have  p  >  1/2-  Tahle  2  gives  the  value  of  B  for  n  =  25,  50,  ® 
for  the  two  Frechet  distributions .  It  also  gives  the  maxinum 
modeling  error  for  the  C  model  vhich  is  an  indication  of 

maximal  excess  modeling  error. 


From  Table  2  we  note  that  the  dependence  structure  exerts  a  large 
effect  on  estimating  the  smaller  of  the  two  component  means  and 
that  either  effect  is  most  exaggerated  fbr  snail  sample  sizes. 

For  K.  1  there  is  very  little  sample  size  effect  on  the  modeling 
error.  For  K  strictly  bigger  than  one  the  maximum  bias  under  the 
Gunbel  C  model  decreases  with  K  and  the  correlation  at  which  this 
maximum  is  attained  also  decreases  to  Q.  Figures  4A-4C  are  plots 
of  B  as  n  ♦  «  for  K  =  3/2,  1,  2/3,  respecitvely.  Figures  5A-5C 
are  plots  of  B  for  n  =  10  for  K  =  3/2,  1,  2/3,  respectively. 


TABU  2 

RELATIVE  R0DELIN6  ERROR  IN  ESTIRATINS  THE  (CAM 


LOVER 

6URSEL  C 

LONER 

SURSEL  C 

K  BOUND 

RH0-1  m  HAS 

BOUND 

RHQ-l  RH0  81A! 

N-INFINITT 


-1*7.44  t««H  1.000  ♦♦♦«♦♦ 
-142.78  •«*♦♦♦  1.000  *♦♦♦♦♦ 
-123.37  ♦♦**«♦  1,000  ****** 
-107.23  ♦«♦»♦*  1.000  «ihm 
-VI.**  ******  1.000  ****** 
-91.71  «♦♦♦♦*  1.000  ****** 
-71. <3  ♦•*♦*♦  1.000  *♦♦♦*♦ 
-*1,78  ******  1.000  ****** 
-32.14  ******  1.000  *«♦♦♦« 


-80.07  ♦♦♦♦♦♦  1.000  •♦«*♦♦ 
-74.13  ******  1.000  ****** 
-74.48  ♦♦♦♦♦♦  1.000  ••♦♦«♦ 
-49.01  ******  1.000  ♦«««♦• 
-43.34  ******  1.000  •♦«♦«♦ 
-41.34  ♦♦*♦*♦’1.000  ****** 
-37.50  *«♦♦♦♦  1.000  •*♦*•♦ 
-31.00  ******  1.000  *♦♦♦«* 
-47.3!  **♦♦♦♦  t.OOO  «*♦♦«♦ 


-38.44  »♦«♦«•  1.000  ♦*•«♦* 
-37.84  ♦♦**«*  1.000  *♦♦♦♦ 
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TABLE  3 

ASTNPOTIC  HAS  QF  THE  PRODUCT  LIMIT  ESTIMATOR 
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Figures  SA-5C  for  p  =  .25,  SA-6C  for  p  =  .5  and  7A-7C  for 
p  =  .75  are  plots  of  &^(p)  for  the  5  models  and  k  =  3/2,  1,  2/3. 
As  in  the  previous  figures  one  can  see  that  for  even  a  small 
departure  from  independence  the  relative  effect  of  dependence 
can  be  quite  large. 


4 . 2  Nonparametric  Estimation 


A  second  approach  to  the  problem  of  estimating  component 
parameters  is  via  the  nonparametric  estimator  of  Kaplan  and 
Meier  Cl  9 58).  Investigators  vho  routinely  use  nonparametric 
techniques  may  take  this  approach  in  hopes  of  obtaining  estimators 
that  are  robust  with  respect  to  the  assumption  of  exponential  itry . 
However  this  estimator  is  not  necessarily  robust  to-  the  assumption 
of  independence. 

The  product  limit  estimator,  assuming  independent  risks  is 

constructed  as  follows.  Suppose  that  n  systems  are  put  on  test 

and  let  rf, ,  — ,  r-  be  the  ranks  of  the  ordered  n.  failures 
jll  in^  i 

from  cause  i,  -•••  )»  among  all  n  order  lifetimes. 

The  estimator  of  the  ocmponent  reliability  for  the  i**1  component  is 


i 


-.2.1.  x<  xi(l)  • 

j(i.x)  x  >  x(i(1) 


j=l  n  -  r^  ♦  1 

Uiere  j(i,x)  is  the  largest  value  of  j  for  thich  x^jj  4  x.  This 
estimator  is  asympotically  unbiased  when  the  component  lifetimes 
are  independent. 

When  the  risks  are  dependent  Klein  and  Moeschberger  (1983) 
show  that  S^(t)  is  not  estimating  the  marginal  component  relia¬ 
bility,  but  rather  it  is  estimating  consistently  another  survival 
function 


4.2.2 


ff^Gc) 


d  QiCt) 
Tit) 


) 


there 


V 


F(t)  =  P (minimum  C^,^)  >  t)  and  QiCt)  = 

PCminGC^)  <  t,  minO^.Xj)  ’  X.),  i  *  1,  2.  Expressions  for 


K-Ct)  for  the  five  models  of  interest  -ire  given  in  Appendix  3. 

A  measure  of  the  affect  of  dependence  in  using  the  Product 
limit  estimator  with  dependent  risks  is  &^Cp)  =  C^Ct^J  -  p)/p 
where  t^  is  the  time  vhere  the  true  component  reliability  is  p. 
A^Cp)  is  again  only  a  function  of  k  =  ?or  t^ie  uPPer 

Frechet  distribution 

dLCp)  =  p_1  -  1  for  k  <  1 

-  1  for  k  =  1  since 

Q  for  k  >  1 

for  k  <  1  H^(t)  =  1  for  all  t  since  the  first  component  never 
fails,  while  for  k  >  1  all  failures  are  due  to  the  first  component. 
For- those  models  with  correlation  spanning  the  range  [0  -  1],  A^(p) 
is  increasing  for  correlations  less  than  p*  and  decreasing  for 
correlations  greater  than  p*  viien  k  >  1.  For  the  lower  Frechet 

distribution  H^Ct^)  =  exp  -  f1  JL _  du  for  p  >_  Cl-Y) 

'  p  u  ♦  u  -  1 

Q  otherwise. 

Table  3  shows  the  value  of  A^(p)  x  100%  for  p  =  .7,.  .5,  .3  for  the 
two  Frechet  distributions.  For  k  >  1,  the  maximum  value  under 
the  Gunbel  C  model  is  also  given.  .As  in  the  parametric  estimation 
problem  the  largest  errors  are  incurred  when  k  <  1.  In  all  cases 
the  effect  of  a  departure  from  independence  is  the  largest  when 
p  is  snail  (i.e.  for  large  t) .  The  effect  decreases  as  k 
increases  reflecting  the  fact  that  vhen  k  X  ^  >  >  X  ^  the  majority 
of  the  system  failures  are  due  to  the  failure  of  the  first 
component. 


5.  Conclusions 

The  results  presented  in  this  paper  show  that  for  all  five 
bivariate  exponential  models  one  nay  be  appreciably  misled  by 
falsely  as suning  independence  of  component  lifetimes  in  a  series 
system.  The  amount  of  error  incurred  in  modeling  system  reliability 
not  only  depends  upon  the  correlation  between  component  lifetimes 
but  also  on  the  level  of  system  reliability.  The  error  in 
modeling  nean  system  life  similarly  depends  upon  the  correlation 
and  the  length  of  mean  system  life.  Both  quantities  depend  on  the 
relative  magnitudes  of  the  parameters. 

For  the  dual  problem  of  estimating  component  reliability 
based  on  data  from  a  series  system,  it  appears  that  departures 
from  independence  are  of  greater  consequence.  Both  parametric 
and  nonparametric  estimators  of  relevant  component  parameters 
are  inconsistent.  Bias  increases  dramatically  as  the  correlation 
gets  further  from  zero.  However,  the  five  models  do  not  exhibit 
appreciable  differences  in  bias  and  mean  squared  error  as 
correlation  changes.  This  suggests  that  these  models  may  belong 
to  a  large  class  of  bivariate  exponential  distributions  which 
possesses  the  properties  exhibited  here. 

ACKNOWLEDGEMENT 

Research  sponsored  by  the  Air  Force  Office  of  Scientific 
Research  under  contract  AFOSR-82-0307 . 


REFERBJCES 


Barlow,  R.E.  and  Proschan,  ?,  0.975).  Statistical  TEtedry  of 
Reliability  and  Life  Testing.  Holt  Rinehart ,  and  Winston. 

Ba.su ,  A.P.  0.981).  Edentifiahility  problems  in  the  theory  of 
competing  and  complementary  risks  -  a  survey.  Statistical 
Distribution  in  Scientific  Vfork  (Taillie,  Patil,  and  Baldesaari , 
Eds.),  Dorerecht,  Holland:  Reidel  Publishing  Co.,  335-34 8 . 

Basu,  A.P.,  and  Klein,  J.P.  (1982).  Some  Recent  Results  in 

Competing  Risks  Theory.  Survival  Analysis,  Crowley  and  Johnson, 
eds.  Institute  of  Mathematical  Statistics,  Hayward,  CA,  216- 
229. 

Bimbaum,  2.W.  (1979).  On  the  Mathematics  of  Competing  Risks. 

U.S.  Dept,  of  HEW.  Publication  No.  79-1351. 

Boardman,  T.J.  and  Kendell,  B.J.  (1970).  Estimation  in  compound 
exponential  failure  models.  Tachnometrics  12,  891-900. 

Clayton,  D.G.  (1978).  A  Model  for  Association  in  Bivariate  Life- 
tables  and  its  Application  in  Epidemiological  Studies  of 
Familial  Tendency  in  Chronic  Disease  Incidence,  Biometrika 
65,  141-151. 

David,  H.A.  and  Moeschberger ,  M.L.  (1978).  The  Theory  of 
Competing  Risks.  Griffin,  London. 

Downton,  F.  (1970).  Bivariate  exponential  distributions  in 
reliability  theory.  J.  Roy.  Statist.  Soc.  B.  32:  408-417. 

Easterling,  R.G.  (1980).  Book  review  of  Z.W.  Bimhaum  monograph. 
Tachnometrics  22,  131-132. 

Frechet,  M.  Sur  les  Tableaux  de  Correlation  Dont  les  Marges  Sont 
Donnees.  Annales  de  1 ' Uni vers it e  de  Lyon.  Section  A,  Series 
3,  14,  53-7 T. 

Gumbel,  E.J.  (196Q) .  Bivariate  exponential  distributions'. 

J.  Aroer.  Statist.  Assoc.  55,  S98-7Q7. 

Hutchinson,  T.P.  (1981).  Compound  Ganna  Bivariate  Distributions. 
Metrika,  28,  263-271. 

Kaplan,  E.L.  and  Meier,  P.  (1958).  Non  parametric  estimation  from 
incomplete  observations.  J.  Amer.  Statist.  Assoc.  53,  457-481. 


uf  I.V.  Atk 


*  rw  v  «  ■  w  in  i 


Klein,  J.P.  and  Moeschberger,  M.L.  C1984).  Asymptotic  bias 

of  the  product  limit  estimator  under  dependent  competing  risks. 
IAPQR  Transactions,  2,  1-7. 

Lagakos,  S.W.  0979).  General  rigjvt  censoring  and  its  impact  on 
the  analysis  of  survival  data.  Biometrics  35,  139-156. 

Lindley,  D.  and  Singpurualla,  N.  0984).  Multivariate  distributions 
for  the  reliability  of  a  system  of  components  sharing  a  cotrmon 
environment.  Presented  at  the  International  Conference  on 
Reliability  and  Quality  Control. 

Mann,  N.R.  and  Grubbs,  F.E.  (1974).  Approximately  optimum  con¬ 
fidence  bounds  for  system  reliability  based  on  component  test 
data.  Technometrics  16,  335-347. 

Mann,  N.R. ,  Schafer,  R.E.,  and  Singpurwalla,  N.D.  (1974).  Methods 
for  the  Statistical  Analysis  of  Reliability  and  Life  Data. 

Wiley,  New  York. 

Miyaraura,  T.  (1982)  ■  Estimating  component  failure  rates  from 
combined  component  and  systems  data:  exponentially  distri¬ 
buted  component  lifetimes .  Technome tries  24:  313-318. 

Miller,  D.R.  (1977).  A  note  on  independence  of  multivariate  life¬ 
times  in  competing  risks.  Ann.  Statist.  5:  576-579. 

Moeschberger,  M.L.  and  David,  H.A.  (1971).  Life  tests  under 

competing  causes  of  failure  and  the  theory  of  competing  risks. 
Biometrika  27,  909-933. 

Moeschberger,  M.L.  and  Klein,  J.P.  (198*0.  Consequences  of 
departures  from  independence  in  exponential  series  systems. 
Technometrics  26:  277-284. 

Oakes,  D.  (1982).  A  model  for  association  in  bivariate  survival 
data.  J.  Boy.  Statist.  Soc .  B.  44:  414-422. 

Peterson,  A.V.,  Jr.  (1976).  Bounds  for  a  joint  distribution 

function  with  fixed  sub-distribution  functions:  Application  to 
competing  risks.  Proc.  Nat.  Acad.  Sci.  73:  11-13. 

Williams,  J.S.- and  Lagakos,  S.W.  (1977).  Models  for  censored 
survival  analysis:  Constnat-sum  and  variable-sum  models. 
Biometrika  64:  215-224. 


civ  o.'i 


Oakes 


1 


i 


>»v5 


-  Formulas  for  p  =  P(X^  <  X^) 


Gumbel  A  - 


PCXj^  <  X?)  =  1/2 


*  ^f)*  (£$■  »X2  > 


where  $  C  • )  is  the  survival  function  of  a  standard  normal 


random  variable. 


(A2.1) 


Gumbel  3  -  P(XL  <  X?)  =  xL 


♦  <+0  A‘2(A1“X2) 


CX1*X2) 


^XyUj^  2X2)UX1+X2) 


Gumbel  C  - 


P(XX  <  x2) 


(x™n™) 


(A2.2) 


(A2.3) 


Downton  -  PCX.  <  X  )  =  2x,X.(l-p) 


(/(x1n2)  -  «*p  x1x2)(x1-x2  ♦  -  4oxLx2) 


(A2.4) 
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Found  numerically  due  to  no  close  fo^n  solution  for 
F(t). 


H^x)  s  exp 
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Asymptotic  Modeling  Error  in  Estimating  for  K  =  1 
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FIGURE  SA 


Small  Sample  Size  (N=10)  Modeling  Error  in  Estimating  for  K  =  l.S 
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FIGURE  6A 

Asymptotic  Error  in  the  Nonparametric  Estimator-  of  the 
First  Component  Survival  Function  at 

t  =  .25,  K  =  1.5. 
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Asymptotic  Error  in  the  Nonparametric  Estimator  of  the 
First  Component  Survival  Function  at 
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ABSTRACT 

A  random  environmental  effects  model  is  proposed  for  competing  risks  experiments.  The 
model  assumes  a  random  stress,  Z,  which  changes  the  scale  parameter  of  each  of  the  assumed 
Weibull  times  to  occurrence  of  the  risks.  Some  general  properties  of  the  model  are  discussed,  and 
specific  properties  for  a  Uniform  or  Gamma  stress  model  are  presented.  Estimation  of  parameters 
under  the  Gamma  stress  model  is  considered,  and  a  new  estimator  based  on  the  scaled  total  time  on 
test  transform  is  presented. 

INTRODUCTION 

The  problem  of  competing  risks- arises  naturally  in  a  number  of  engineering  or  biological 
experiments.  In  such  experiments,  for  some  items  put  on  test,  the  primary  event  of  interests  (such 
as  death,  component  failure,  etc.)  is  not  observable  due  to  the  occurrence  of  some  competing  risk 
of  removal  from  the  study  (such  as  censoring,  failure  from  a  different  component,  etc.). 

Competing  risks  arise  in  an  engineering  context  in  analyzing  data  from 

(a)  series  systems, 

(b)  field  tests  of  equipment  with  a  fixed  test  time  and  a  random  or  staggered  entry  into 
the  study,  or 

(c)  systems  with  multiple  failure  modes. 

Competing  risks  arise  in  biological  applications  in  analyzing  data  from 


(a)  clinical  trials  with  a  fixed  trial  duration  and  staggered  entry 

(b)  clinical  trials  with  some  patients  withdrawing  from  the  trial  prior  to  response 

(c)  studies  of  the  time  to  death  from  a  variety  of  causes 

A  common  assumption  made  in  analyzing  competing  risks  experiments  is  that  the  potential 
(unobservable)  times  to  occurrence  of  the  competing  risks  are  independent  This  assumption  is  not 
testable  due  to  the  identifiability  problem.  That  is,  for  any  dependent  competing  risks  model,  there 
exists  an  independent  competing  risks  model  which  yields  the  same  observables.  (See  Basu  and 
Klein  (1982)  for  details.)  However,  Moeschberger  and  Klein  (1984)  show  that  an  investigator  can 
be  appreciably  misled  in  modeling  competing  risks  by  erroneously  assuming  indepencence. 

In  this  paper  we  present  a  model  for  dependence  between  the  various  risks  by  assuming  that 
dependence  is  due  to  some  common  environmental  factor  which  effects  the  potential  times  to 
occurrences  of  each  risk.  In  section  2  we  present  the  model  and  study  its  properties  for  bivariate 
series  and  parallel  systems.  In  section  3,  we  consider  estimation  of  the  model  parameters  for 
competing  risks  systems. 

2.  THE  MODEL 

For  simplicity  we  shall  consider  the  problem  of  bivariate  systems  and  discuss  our  model  in 
terms  of  engineering  applications.  We  assume  that  under  ideal,  controlled  conditions,  as  one  may 
encounter  in  the  laboratory  in  the  testing  or  design  stage  of  development,  the  time  to  failure  of  the 
two  components,  to  be  linked  in  a  system,  are  Xq  and  Yq.  We  suppose  that  under  these 

conditions,  Xq,Yq  have  survival  functions  Fq,  Gq  on  [0,  <»).  We  assume  that  both  Xq  andYg 

follow  a  Weibull  form  with  parameters  (T|p  Xp  and  (  t^,  X^),  respectively.  That  is,  Fq(x)  =  exp(- 

Xjx^l).  The  Weibull  distribution,  which  may  have  increasing  ( T|  >1),  decreasing  (T|  <  1)  or 

constant  failure  rate  (T|  =  1)  has  been  shown  experimentally  to  provide  a  reasonable  fit  to  many 
different  types  of  survival  data.  (See  Bain  (1978)).  We  now  link  the  two  components  into  a 
system  in  such  a  way  that  under  ideal  lab  conditions  the  two  components  are  independent 

Now  suppose  that  the  above  system  (Xq,  Yq)  is  put  into  operation  under  usage  conditions. 

We  suppose  that  under  such  conditions  the  effect  of  the  environment  is  to  degrade  or  improve  each 
component  by  the  same  random  amount  That  is,  the  effect  of  the  environment  is  to  select  a 
random  factor,  Z,  from  some  distribution,  H,  which  changes  the  maginal  survival  functions  of  the 
two  components  to  Fq'-  and  Gq  .  A  value  of  Z  less  than  one  means  that  component  reliabilities  are 
simultaneously  improved,  while  a  value  of  Z  greater  than  one  implies  a  joint  degradation.  The 
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resulting  joint  reliability  of  the  two  components'  lifetimes,  (X,Y)  in  the  operating  environment  is 

F(x,y)  =  E[exp(-Z(X1xTU+  X2yT,2  )].  (2.1) 

This  model  has  been  proposed  by  Lindley  and  Singpurwalla  (1984)  in  the  reliability  context 
when  Fq,  Gq  are  exponential  and  H(  )  follows  a  gamma  distribution.  This  basic  dependence 

structure  was  also  proposed  by  Clayton  (1978)  to  model  associations  in  bivariate  survival  data,  and 
later  by  Oakes  (1982)  to  model  bivariate  survival  data.  Hutchinson  (1982)  proposed  a  similar 

model  when  H(  )  has  a  gamma  distribution  and  Fo(t)  =  Go(t)  =  exp(-t^). 

The  model  described  above  for  a  general  distribution  of  the  environmental  stress  has  a 
particular  dependence  structure  which  we  summarize  in  the  following  lemmas. 

Lemma  1.  Let  (X,Y)  follow  the  model  (2.1)  where  Z  is  a  positive  random  variable  with  finite 

r  s 

( - + - )^  inverse  moment  Then 

^1 

*r/TU  -s/rj2  -(rim  s/m) 

E(XrYs)  =  X1  X2  HI  +  r/rij)  T(1  +  s/tl2)  E(Z  )  (2.2) 

The  proof  follows  by  noting  that,  given  Z  =  z,  (X,Y)  are  independent  Weibulls  with  parameters 

-rAi!  -r/ri! 

(Tlj,  z)  and  (t^,  z),  respectively  and  E(Xr|Z=z)  ==  Xj  z  T(l+  rirti)  with  a  similar 

expression  for  Ys.  When  the  appropriate  moments  exist,  we  have 

1  bt\ 

(A)  E(X)  =  E(Xq)  E(Z  S, 

-1/rii  -lAl!  , 

(B)  V(X)  =  E(Xo2)Var(Z  )  +  E(Z  )2  VarfXo), 

-1/Tlj  -1/T i, 

(C)  Cov  (X,Y)  =  E(Xq)  E(Yq)  Cov(Z  ,  Z  )  which  is  greater  than  0. 

If  ill  =  r|2  =  t|  then  the  correlation  between  (X,Y)  is 

r(l+  1/n)2  Var(Z-1/TI) 

P= - 

VarOZ*1^)  r(l+2ril)+  (T(l+2/n)  -  T(1+1/ti)2)  E(Z’l71)2 

In  this  case  the  correlation  is  bounded  above  by  r(l+l/rj)2  /  r(l+2/T|).  Figure  1  shows  the 

maximal  correlation  as  a  functeion  of  r|  for  ti  e  (0, 10).  Note  that  this  maximal  correlation  is  an 
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increasing  function  of  rj.  One  can  also  show  that  F(x,  y)  is  positive  quadrant  dependent  for  any 


^  1’  t'2* 

Exact  expressions  for  competing  risks  quantities  of  interest  can  be  computed  when  a 
particular  model  is  assumed  for  the  distribution  of  Z.  We  shall  consider  the  gamma  and  uniform 

models.  Consider  first  the  gamma  model  with  hz(z)  =  ba  za' 1  exp(-bz)/T(a),  z  >  0.  For  this 
model,  the  joint  survival  function  is 


ba 


F(x,y)  = 


(2.3) 


[b+ 


r\i  ill, 

Xix  +  T^y  la 


which  is  a  bivariate  Burr  Distribution  (see  Takahasi  (1965)),  the  marginal  distributions  are 
univariate  Burr  distributions  with 


*1/TU 

E(X)  =  (Xt/b)  ni+l/Hi^a-  lAli)  Ha),  if  a  >  1/q, 


-2ai1  ro+^^rca-^!)  ra+i/n^na-i/tij) 

Var (X)  =  (  Xi/b)  { - -  [ - ]2},  if  a  >  2/1^ 

r  (a)  r(a) 


with  similar  expressions  for  E(Y),  Var(Y).  The  covariance  of  (X,Y)  is 


-1/Tji  -1/T|2  na-l/Tjt-l/ri^)  r(a-l/ti2)  Ha-l/TH) 

cov(x,Y)  =  (Xj^/b)  (k2/ b)  ra+i/ndni+i/i^K - - - ■> 


r(a) 


r(a) 


for  a  >  l/r|  j  +  1/T12-  For  the  gamma  model,  the  reliability  function  for  a  bivariate  series  system  is 
given  by 


^1  ^2  a 
Rj(t)  =  (l+fl^/bjt  +  (X2/b)t  )*a, 


(2.4) 


and  for  a  parallel  system  by 


Til 


rl2. 


^2 


T12 


Rp(t)  =  (l+(X1/b)t  +  (1+  (X2/b)t  )*a-  (l+(  Xj/bJt  +  (X2/b)t  )‘a  (2.5) 

Figures  2A-E  are  plots  of  the  series  system  reliability  for  X^=  1,  ^  =  2  and  several  combinations 
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of  Tip  ri2*  Each  figure  shows  the  reliablity  for  a  =  1/2, 1, 2, 4,  and  the  independent  Weibull 

model.  In  all  cases,  b  =  1.  For  these  figures  we  note  that  for  fixed  Xp  X2  "Hi5  ^2*  **  the  series 

system  reliability  is  a  decreasing  function  of  the  shape  parameter  a.  Figures  3A-E  are  plots  of  the 
parallel  system  reliability  (2.5)  for  the  above  parameters.  Again,  the  reliability  is  a  decreasing 
function  of  a.  Also  in  both  the  series  and  parallel  system  reliability,  the  shape  of  the  reliability 
function  is  quite  different  from  that  encountered  under  independence. 

The  gamma  model  is  a  reasonable  model  for  the  environmental  stress  due  to  its  flexibility  and 
the  tractability  of  the  model  in  obtaining  close  form  solutions  for  the  relevant  quantities  and  in 
estimating  parameters.  However,  in  some  cases,  such  as  when  the  operating  environment  is 
always  more  severe  than  the  laboratory  environment,  the  support  of  H  may  be  restricted  to  some 
fixed  interval.  A  possible  model  for  such  an  environmental  stress  is  the  uniform  distribution  over 
[a,b].  For  this  model,  the  joint  survival  function  is 


Til  ^2  ^1  ^2 

[exp  (-b(Xj  x  +  X2y  ))  -  exp(-a(X1x  +  X^y  ))] 
F(x,y)  =  — - - 

Hi  112 
(b-aXXjx  +  X2y  ) 

-l/ill  Olrl)/Tii  OlrD/m 

E(X)=Xj  ra+l/TliJTija,  -a  )/{(Tirl)(b-a)] 

=  JlnCb/aVtX^b-a)] 


(2-6) 

if  ill  ^  1 

if  T|l=l, 


•2/  T| !  ('n1-2)/'n1  (tj1-2)/t|i 

Var(X)  =  q1Xl  {  HI +2/^) rjj  (b  -a  1 


-  (T|i-1/  T|i)  (Tli-l/Tl!). 

-  r  (l-t-l/Ti!)2  ri!  (b  -a  )2  if  rjj  *  1,2 

”  (Tl1-l)2(b-a) 


2/(X^2ab)  -  fln(b/  a)2/[(b-a)  Xj]2 


V1! 


2n(b/  a) 


(bl/2  +  al/2)i 


if  11!=  1 


ifTli=  2 
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FIGURE  2  E 

SERIES  SYSTEM  RELIABILITY  UNDER  GAMMA  (A,  1)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

*1  =  1  *0.  *2=2 . 0 ,  ’>1  =  1/2,  ’>2=2. 0. 
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FIGURE  3  R 

PARALLEL  SYSTEM  RELIABILITY  UNDER  GAMMA  (A,  1)  MODEL 
FOR  THE  EN  r'IRQNMENTAL  STRESS. 

Xi  =  1.0,  *2=2.0,  ^i  =  l.Q.  ’(2=1.0. 
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FIGURE  3  C 

PARALLEL  SYSTEM  RELIABILITY  UNDER  GAMMA  (A, 1)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

>M  =  1.0,  >2=2.0,  *>i=2.  Q.  *2=2.0. 


FIGURE  3  D 

PARALLEL  SYSTEM  RELIABILITY  UNDER  GAMMA  (A,  1)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

>>1  =  1.0,  >>2=2 . 0 ,  *>i  =  l/2.  *12=1/2. 
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*11*12 
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Tli-l  Tlj-1  tl2-l  T12-l 

*ll  *U  *12  *12 

(b  -a  )(b  -a  ) 


r(l+l/ni)r((2Tirl)/tii)  An(b/a) 
1/Hi  Tli^Tli-l)  (b-a) 

Xi  X2 


}if  rit  1,  T12#!.  1/*U  +1/Tl2,tl 


1-1  “Hi-1 

-  m  Tii 

Til2  (b  -a  )  1/Tli  1/Tl! 
— - , - (b  -a  )] 


Oi  rl)  (b-ar 


*li  "Hi*  *  *li  ‘  *li*  *li  *li’  *  *li "  *li' 


if  l/ni+  1/rip  1 


Tli-l  Tli-l 


TliTli- 

r(l+l/Tli)  TljTlj-  (b 

-  [ - 

1/Hi  ( Tjj  Tij.  -  rij  -  TijOCb-a) 
Xi  \ 


*li  Tli 

)  rji  (b  -a  )  Un(b/a) 


(Tlj-l)  (b-a)" 


if  “Hi  **  l.rii—  1 


1  1  lln(b/a)2 

- [ - - - *—  ]  if  Tl!  =ti2=  1 

(Xj  X2)  (ab)  (b-a)2 

For  this  model,  the  reliability  function  for  a  series  system  is 


R$(t)  =  [exp(-b( 


Tli  ri2  ’ll  *12 

X^t  +X2t  )-exp[-a(Xit  +X2t  )]] 


„  Bi  .  t\2 
(b-a)(  Xjt  +  Xj  t  ) 


and  for  a  parallel  system  is 


'll  Hi  ^2  'l2 

Rp(t)  =  [exp(-b(Xjt  )  -exp(-a  Xjt  )  +  [exp(-b  fyt  )  -exp(-a  ) 


^2. 


-Rs(t) 


(b-a)  Xjt 


^1 


(b-a) 


^2 


Figures  4A-E  show  the  reliability  for  a  series  system  and  figures  5A-E  for  a  parallel  system 

under  the  uniform  model  for  various  combinations  of  Xj,  ti  p  tj  2  a>b-  Notice  that  when  A  = 

.25,  B  =  .75,  which  corresponds  to  an  operating  environment  which  is  less  severe  than  the  test 
environment,  the  system  reliability  is  greater  than  that  expected  under  independence,  while  when 
(a,b)  =  (1.25, 1.75)  or  (1.,  2),  which  corresponds  to  an  environment  more  severe  than  the  test 
environment,  the  system  reliability  is  smaller.  Also  when  the  (a,b)  contains  1,  which  corresponds 
to  an  environment  which  incurs  the  possibility  of  no  differential  effect  from  that  found  in  the 
laboratory,  there  is  little  difference  in  the  dependent  and  independent  system  reliability. 


3.  Estimation  of  Parameters  Under  Gamma  Model 


Consider  the  model  (2.3)  with  rj  ^  =  tj2  =  Tl-  For  this  model,  the  reliablity  for  a  series  system 


is 


R$(t)  3  U  + 


(X1+X2) 

b 


t^-a. 


(3.1) 


Notice  that  this  model  depends  only  on  two  parameters  9  =  (Xj+  X^fb  and  a  so  that  if  we  had  data 
only  from  systems  on  test  in  the  operating  environment,  the  only  identifiable  parameters  are  a,  9, 
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FIGURE  4  D 

SERIES  SYSTEM  RELIABILITY  UNDER  UNIFIR, B)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

*1=1.0,  *2=2.0,  *2=1/2. 
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FIGURE  5  fl 

PARALLEL  SYSTEM  RELIABILITY  UNDER  UNIF(A.B)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

Xjsl.Q,  X2  =  2.0,  •’71=1.0.  ^2=1 


FIGURE  5  B 

PARALLEL  SYSTEM  RELIABILITY  UNDER  UNIF(A.B)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

>1  =  1.0,  >2=2.0,  *»i=  l .  0,  n2=2.0. 
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FIGURE  5  C 

PARALLEL  SYSTEM  RELIABILITY  UNDER  UNIF(A.B)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 
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t],  not  Xj,  %2>  "H*  a>  b.  However,  in  many  instances  we  have  extensive  data  on  the  performance  of 


the  components  in  the  lab  under  ideal  operating  conditions  so  that  one  may  consider  Xj,  r\  to  be 

known  based  on  estimates  from  this  data.  We  shall  focus  on  the  problem  of  estimating  0  and  a, 
based  on  data  on  the  system  failure  times  collected  in  the  operating  environment  Let  4, ....  tn  be 

the  failure  times  for  n  such  systems  put  on  test,  and,  let  w^  =  t^ ,  i  =  1, n. 

Prior  to  attempting  to  estimate  (a,  0),  we  would  like  to  check  if  the  model  (3.1)  is  feasible. 

A  graphical  check  of  this  model  can  be  done  through  the  scaled  total  time  on  test  (STTOT)  plot  of 
Barlow  and  Campo  (1975).  The  STTOT  for  W  is 

,F‘(t) 

/  R  s(t)dt 

Gw(t)  - - 7 - =  l-(l-t)(a'1)/a  for  a  >  1.  (3.2) 

fHd 

}  Rs(t)dt 

0 

Note  that  (3.2)  depends  only  on  a.  Figure  6  shows  the  form  of  the  STTOT  for  several  values  of  a. 

Notice  that  for  all  a,  the  STTOT  is  below  the  45°  line  (which  corresponds  to  exponential  system 
life)  since  the  hazard  rate  of  the  series  system  is  decreasing.  Let 
i 

Tn(W(i))  =  IW(j)  +  (n-i)W(i),  (3.3) 

where  ^ ...  £  W^  are  the  ordered  systems  failure  times  be  the  total  time  on  test  at 

The  empirical  STTOT  plot  then  plots  (i/n,  Tn(W^p/T(W(n)))  which  can  be  compared  to 
figure  6  for  a  graphical  check  of  the  model.  Also,  crude  estimates  of  a  can  be  obtained  by 


’  _"*•  _  *  «"*»  V  S  *  -V  J'’  -  ■  •  ’  »  *■  »  ’  k  "  *  *■  »  ■' ’  »  ■*  ta  *  -  *■  r 
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comparing  the  empirical  and  theoretical  STTOT  plots.  When  there  is  no  random  environmental 

effect  and  the  components  are  independent,  then  the  empirical  STTOT  plot  should  look  like  the  45° 

line.  Also  as  a  tends  to  infinity  this  plot  approaches  the  45°  line. 

We  now  consider  several  estimates  of  a  and  0.  The  log  likelihood  for  the  model  (3.1),  based 
on  a  sample  of  size  n,  is 

n 

L(a,0)  =  nHna  +  nlln0-  (a+l)Z  lln  (1+0  W:)  (3.4) 

i=l 


so  that 


n 

9/9a  L(a,0) «  n/a  - Z  lln  (1  +  0  W:) 
i=l 


n 

and  9/90  L(a,0)  =  n/0  -  (a+1)  Z  w:/(l+  0w;) 

i=l 

For  (3.5)  we  note  that  the  maximum  likelihood  estimator  of  a 

n 

amle  ~ - 

n 

S  fln  (1  +  0  Wj) 
i=l 


and  the  maximum  likelihood  estimator  of  0  is  the  solution  to 


n  n  Wj 

-  -( - +1)(I - )  =  0. 

0  I  An  ( 1  +  0  Wj)  l  +  0Wj 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


n 

One  can  show  that  0  is  positive  if  n  Z  w 

i=l 


.2 

l 


n 


>  2(  Z  w:)2. 

i=l 


(3.9) 


In  such  case  0-^-  is  obtained  by  solving 


(3.8)  numerically. 


A  second  estimator  of  (a,  0)  is  the  method  of  moments  (rame).  Since  E(W)  =  [  Q(a- 1)]_1 

\ 

and  E(W2)  =  2[  02(a-l)(a-2)]-1  where  a  >  2,  we  have 


amme=  1  + 


(3.10)  and  0^6  = 


I  wi2-2(Swi)2 
I  Wj(I  Wj)2 


(3.11) 


provided  that  (3.9)  holds.  If  (3.9)  does  not  hold,  then  this  estimator  does  not  exist 

A  third  estimator  was  suggested  by  Berger  (1983)  in  a  different  context  He  suggested 

estimating  0  a  modified  methods  of  moments  estimator  0^  =  (a  w)"1 ,  (3. 12) 
where  w  =  X  w-/n. 


which  is  used  as  the  true  value  of  0  in  the  likelihood  (3.4)  so  that  the  estimator  of  a  is  the  solution 
to 

wt  (a+1)  wj 

-  X  fin  (1+ - )  + - ,1 - 0  (3.13) 

aw  waz  l+wj/(aw) 

A  final  estimator  is  based  on  the  STTOT  plot  Let  Cj  =  fln(l-i/n)  and  Dj  = 
fln(l-Tn(W^p/Tn(W(-np),  i  =  1, ...,  n-1.  If  (3.2)  holds,  then  we  should  have 
Di-(l-l/a)Ci,i=  l,...,n-l,  (3.14) 

so  the  value  of  a  which  minimizes 

n-1  - 

£  (Dj  *  (1-1/a)  Cjr  is  a  reasonable  estimator  of  a 
i=l 

Zq2 

The  resulting  estimator  is  aie  = - * -  (3.15) 


I 

;gj 


"M 


which  is  in  the  parameter  space  if  X  Cj2  >  X  Cj  Dj.  A  better  estimator  should  be  obtained  by 
weighting  the  Dj's  differently  since  for  i  <  j,  Var  (Dj)  <  Var  (Dp.  The  variance  of  Dj  depends  on 
the  unknown  parameter  a  so  we  weight  by  the  variance  of  Dj  computed  under  an  assumed 
exponential  distribution.  The  variance  of  Dj  in  that  case  is 


i  1 

Vi  =  X  - , 

j=l  (n-j> 


i  =  I,-,  n-1 


(3.16) 


so  that  the  weighted  least  squares  estimator  of  a  is 


^Is  “■ 


C,D, 
2— > 
vi 


ifXCj2/Vj>XCjDj/Vj. 


(3.17) 


Once  we  have  obtained  a  by  either  of  the  two  least  squares  estimators,  we  substitute  this  value  into 
(3.6)  and  solve  this  equation  numerically  for  8^  or  0wjs. 

The  condition  X  Cj  2/V;>Iq  Dj/Vi  includes  a  few  more  possible  samples  than  the 
condition  (3.9)  for  the  other  three  estimators.  However,  those  samples  which  satisfy  X  Cj2/Vj  > 

X  Ci  Dj  /Vj  for  which  (3.9)  fails  to  obtain  yield  very  large  estimates  of  9.  Since  a  reasonable 
model  for  T  when  9  and  a  are  not  estimable  is  the  independent  Weibull  series  system  which  has 
system  reliability  very  close  to  (3. 1)  when  a  is  very  large,  this  is  not  a  problem.  Figures  7a  and  7b 
are  scaled  total  time  on  test  plots  from  two  simuilated  samples  of  size  30  from  (3.1)  with  a  =  3,  0  = 
1.  Looking  at  figure  7a,  we  see  that  the  estimated  scaled  total  on  test  doesn't  look  too  different 
from  the  45°  so  that  an  exponential  model  might  not  be  unreasonable.  For  this  data  set  only  the 
weighted  least  squares  estimator  exists  and  it  yields  awT.s  =  45.33  and  8  =  .0567.  For  the  data  in 
figure  7b  all  estimates  exist,  and  we  have 


SCALED  TOTAL  TIME 

0.40  0.60  0.80  1.00 
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®mle  ~  -93  amle  “  ^.98 

=  amme  =  ^-86 

Qber  -  -720  =  7.02 


01s  =  .739  als  =  3.58 

®wls  =  -^70  ay/ls  =  2-89 

To  study  the  properties  of  these  estimators,  a  small  scale  Monte  Carlo  study  was  performed. 
Random  samples  of  size  n  =  15,  30,  50,  75,  or  100  were  generated  with  +  ^2  =  3,  b  =  3,  so  0 
-  1  and  a  =  2, 3, 5.  1000  samples  were  generated  for  each  combination  of  n  and  a.  The  bias, 
standard  deviation  of  the  estimates  and  n,  the  number  of  samples  where  the  estimator  exists  is 
reported  in  table  1  for  a,  table  2  for  0,  and  in  table  3  for  an  estimator  of  the  system  reliability 
obtained  from  (3. 1)  at  tg  =  9.085.  The  true  system  reliability  at  to  is  .8255  when  a  =  2,  .75  when  a 
=  3,  and  .619  when  a  =  5.  Also  reported  in  each  table  is  the  bias  and  standard  deviation  of  the  least 
square  and  weighted  least  square  estimators  when  they  are  restricted  to  those  samples  where  the 
other  estimators  exist 

From  these  tables  we  note  that  Berger's  modified  estimator  performs  very  poorly.  Also  the 
weighted  least  squares  estimator  allows  for  estimation  of  parameters  in  many  more  samples  when  n 
is  small.  In  general  the  maximum  likelihood  estimator  outperforms  the  other  estimators,  however, 
when  the  weighted  least  squares  estimator  is  restricted  to  those  samples  where  the  maximum 
likelihood  estimator  exists,  this  estimator  performs  much  better  when  n  is  small.  The  somewhat 
better  performance  of  the  MLE  in  terms  of  bias  is  deceptive  since  some  of  the  estimates  of  a  are  less 
than  one,  which  implies  that  the  mean  system  reiliability  is  infinite.  Also  the  weighted  least  squares 
estimator  of  system  reliability  seems  to  outperform  the  other  estimators  of  the  system  reliability  in 
spite  of  its  relatively  poor  performance  as  an  estimator  of  0.  Our  recommendation  is  to  use  the 
weighted  least  squares  estimator  since  it  more  often  provides  estimators  of  the  relevant  parameters 
and  is  somewhat  easier  to  compute. 
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1.  Introduction 

The  problem  of  competing  risks  arises  naturally  in  a  number  of  contexts, 
namely  the  modeling  of  series  systems  in  reliability,  the  problem  of  estimation 
with  censored  data  and  the  analysis  of  physical  or  biological  systems  with 
mulciple  failure  modes.  A  common,  uncestable  assumption  is  usually  made  thac 
che  pocendal  failure  times  for  each  risk  are  statistically  independent  (see 
Basu  and  Klein  (1982)).  Moeschberger  and  Klein  (1984)  show  that  an  investiga¬ 
tion  may  be  appreciably  misled  in  modeling  series  systems  reliability  and  in 
escimacing  component  parameters  by  Incorrectly  assuming  independent  component 
lifetimes.  In  this  paper  we  model  dependence  between  components  through  a 
common  environmental  effect  on  each  component.  Such  a  dependence  structure 
has  also  been  suggested  by  Oakes  (1982),  Llndleyand  Singpurvalla  (1985),  and 
Hutchinson  (1981). 

2.  The  Model 

Consider  a  two-component  system.  Suppose  that  under  ideal,  controlled 
conditions,  as  encountered  in  the  testing  stage,  the  times  to  failure  of  the 
two  components  are  XQ  and  Y  .  Under  these  conditions,  XQ,Y_  are  independent 
with  marginal  survival  functions,  FQ  and  Gq.  Now  suppose  cne  two  components 
are  linked  into  a  system  and  exposed  to  the  environment.  The  effect  of  the 
environment  is  to  select  a  random  factor,  Z,  from  a  distribution,  H(z)  ,  which 

Z  Z 

changes  che  marginal  survival  functions  of  the  two  components  to  , 

respectively.  A  value  of  Z  less  than  one  means  chat  component  reliability  is 

improved,  while  a  value  greater  chan  one  implies  a  joint  degradation.  In  che 

sequel  we  assume  chat  X_  and  Y„  follow  a  Weibull  distribution  with  parameters 

w  U  cl. 

( a  '  X  >  ’  and  Thac  is,  ^(x)  “  exp("^xx  *j .  The  resulting  joint 

reliability  of  che  two  components'  lifetimes,  (X,Y),  in  the  operating  environ- 

°X  °Y 

menc  is  F(x,y)  ■  E(exp(-Z(X^x  ♦  X^y  J).  F(x,y)  is  positive  quadrant  dependent. 
Also  E(X)  -  E(X0JE(Z~L/'*)  ;  V(X)  -  E(XQ2)  E(Z"2/Clx)  -(E  (XgjEU"1^0*)) 2  ;  and 

Cov(X.Y)  -  E(X0)E(Yq)  Cov(Z~1/a*,  Z-1^).  The  correlation  is  always  positive 

and  is  bounded  above  by  T(1  l/a)2/r(l  ♦  2/a)  when  ax  •  a  •  a.  Explicit, 

chough  lengthy,  formula  for  che  momenta  of  (X,Y),  system  and  component 
reliability,  and  system  components  can  be  obcained  when  Z  is  assumed  to  be 
elcher  a  uniform  or  gammar  random  variables.  The  gamma  case  leads  to 
bivariate  Burr  distributions. 

3.  Estimation 

The  esc  lame  ion  of  model  parameters  is  carried  out  under  Che  assumption 

that  7.  has  a  gamma  distribution  with  density  hCsJa**”1  exp(-bs).  Since  che 
parameters  X_.  X„  are  not  identifiable  when  only  data  from  either  a  series  or 


parallel  syaca«  la  available,  we  Incorporate  sample  information  obtained 
Independently  on  each  component  under  teat  conditions.  Maximum  likelihood 
and  method  of  moments  estimators  are  obtained  and  their  properties  are  studied 
by  Monte-Carlo  methods  since  no  closed  form  maximum  likelihood  estimates  are 
available. 
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Resume 

Un  modele  pour  les  systemes  dependants  dans  l'analyse  de  la  fiabilite  esc 
examine.  Le  modele  suppose  qua,  sous  les  conditions  ideales,  les  temps  de 
survle  dcs  constituents  du  systems  ont  des  distributions  Uelbull  independences . 
Sous  des  conditions  d'operation  un  facteur  exterieur  aleacoire  affecce  cheque 
constituent  slmultandment  an  multipliant  son  caux  de  hasard  par  une  quant  ltd 
alfatoire.  Les  proprlctes  de  ce  modele  et  l'astlmation  des  parametres  du 
modkle  sont  conslddrds,  k  partir  des  exemples  concrets  du  laboratoirc  et  de 
la  pratique. 


